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OBJECTIVES: 


This  3-year  research  program  focuses  on  micromechanics  and  fracture  mechanics 
analyses  of  cracking  in  functionally  graded  coating/substrate  systems  used  in  turbine  engines 
and  wear-related  applications.  The  objectives  of  the  study  include:  quantifying  the  influence  of 
coating  gradation,  evaluating  the  effect  of  metal  plasticity  and  crack  bridging,  and  predicting  the 
fracture  driving  forces  and  fracture  resistance  curves.  Emphases  is  placed  on  multiple  cracking 
and  delamination  fracture  in  coating/substrate  systems.  The  fundamental  understanding  gained 
from  this  study  will  enable  one  to  select  coating  gradation  at  the  component  design  stage. 


SUMMARY  OF  EFFORT: 

In  the  first  six  months  (Feb.  1995  -  July  1995),  research  was  conducted  to  develop 
micromechanics  models  for  predicting  the  thermomechanical  properties  of  functionally  graded 
coatings.  Coatings  of  the  inclusion/matrix  type  were  considered,  with  the  coating  gradation 
characterized  by  a  power-law  type  relation.  The  elastic  moduli  and  thermal  expansion 
coefficient  of  the  coating  were  predicted  in  terms  of  the  local  volume  fractions  and 
thermomechanical  properties  of  the  phases  in  the  coating.  Distribution  of  thermal  residual 
stresses  in  the  coating/substrate  systems  were  calculated.  A  linear  fracture  mechanics  analysis 
was  carried  out  to  calculate  the  crack  driving  forces  (energy  release  rates  or  stress  intensity 
factors)  for  multiple  cracks  in  the  coating  under  both  thermal  and  mechanical  loads. 

During  the  12-month  period  of  August,  1995  -  July  1996,  a  linear  fracture  mechanics 
analysis  was  carried  out  to  calculate  the  crack  driving  forces  (energy  release  rates)  for 
delamination  cracks  in  the  coating  or  at  the  coating/substrate  interface.  Multilayered  coatings 
and  coatings  of  the  inclusion/matrix  type  were  considered,  with  the  coating  gradation 
characterized  by  a  power-law  type  relation.  A  steady-state  heat  transfer  analysis  was  also 
performed  for  the  FGM  coating/substrate  system  to  calculate  the  increase  of  the  coating 
thickness  owing  to  the  functional  gradation.  A  crack  bridging  model  was  developed  for  the 
crack  tip  shielding  effect  due  to  the  non-uniform  distribution  of  metal  ligaments  in  the  crack 
wake.  The  reduction  in  the  fracture  driving  force  owing  to  crack  bridging  was  quantified. 

From  August  1996  to  July  1997,  research  was  conducted  to  develop  computational 
models  for  bridged  edge  and  internal  through  cracks  in  functionally  graded  coating/substrate 
systems.  A  singular  integral  equation  approach  was  taken  to  solve  the  boundary  value  problem 
in  which  both  the  coating  and  the  substrate  are  assumed  to  be  linear  elastic.  The  crack  driving 
forces  of  edge  and  internal  cracks  in  a  FGM  panel  and  in  a  FGM  coating/metal  substrate  system 
were  calculated  and  the  effect  of  the  coating  gradation  quantified.  An  analysis  was  also  carried 
out  for  cracking  in  TBC/substrate  systems  consist  of  TBC,  FGM  bond  coat,  and  metal 
substrate. 

During  the  period  of  August  1997  -  January  1998,  an  effort  was  made  to  predict  the 
fatigue  behavior  of  advanced  ceramics  and  ceramic  coatings.  A  model  was  developed  for 
fatigue  growth  of  an  edge  crack  in  nontransforming  ceramics  based  on  a  crack  bridging 
analysis.  Systematic  predictions  were  made  of  the  effects  of  initial  crack  length,  load 
frequency,  and  the  bridging  characteristics.  This  work  could  be  further  extended  to  include  the 
effect  of  nonuniform  distribution  of  bridging  grains. 


BRIEF  OUTLINE  OF  RESEARCH  FINDINGS: 

Functionally  graded  coatings  such  as  thermal  barrier  coatings  on  advanced  metals  and 
alloys  can  sustain  multiple  cracking  upon  thermal  and  mechanical  loading.  A  special  feature 


pertaining  to  the  fracture  behavior  of  coating  is  the  coating  gradation  which  can  be  characterized 
by  the  local  volume  fraction  of  metal.  To  guide  the  design  of  coating  gradation  so  that  cracking 
damage  can  be  minimized,  a  fracture  mechanics  study  is  made  of  the  crack  driving  force  for 
multiple  cracks  in  functionally  graded  ceramic/metal  coatings.  The  metal  substrate  and  the 
ceramic/metal  coating  are  taken  as  linear  elastic,  with  the  elastic  properties  of  the  coating 
varying  through  the  film  thickness.  Systematic  finite  element  calculations  are  made  for  the 
energy  release  rate  of  the  cracks  in  the  coating  as  determined  by  the  coating  gradation,  crack 
length,  and  the  crack  density;  both  mechanical  and  thermal  loads  are  considered.  It  is  found 
that compared  with  the  pure  ceramic  coating,  gradation  of  the  coating  can  significantly  reduce 
the  crack  driving  force.  It  is  also  found  that  under  mechanical  loading  the  effect  of  different 
gradations  on  the  crack  driving  force  is  relatively  small.  However,  under  thermal  loading  the 
influence  of  coating  gradation  can  be  significant. 

Delamination  cracking  is  a  major  damage  mode  in  multilayered  composites  including 
FGMs.  To  quantify  such  damage,  a  micromechanics  analysis  is  carried  out  of  delamination 
cracking  in  functionally  graded  coating/substrate  systems.  The  FGM  coating  is  taken  to  be  a 
ceramic/metal  composite  with  its  gradation  characterized  by  the  local  volume  fraction  of  metal. 
The  elastic  properties  thus  change  with  position  in  the  coating  along  the  thickness  direction. 
Linear  fracture  mechanics  analyses  are  made  for  the  strain  energy  release  rate  of  the 
delamination  crack  as  determined  by  the  coating  gradation,  the  thickness  ratio  of  the 
coating/substrate  system,  and  the  elastic  properties  of  the  ceramic  and  metal  phase.  Both  the 
total  energy  release  rate  and  the  stress  intensity  factors  of  the  system  are  predicted.  The  mode 
mixity  can  be  quantified  using  the  analytical  formula  developed.  Buckling  of  the  FGM  coating 
due  to  the  compressive  residual  stresses  is  analyzed  based  on  the  thin  plate  theory;  the  buckle- 
driven  delamination  is  modeled  accordingly.  It  is  found  that  functional  gradation  of  the  coating 
has  a  strong  influence  on  the  residual  stress  distribution,  which  in  turn  changes  the  buckling 
behavior  of  the  coating.  By  using  a  FGM  coating,  the  fracture  driving  force  for  the 
delamination  crack  can  be  significantly  reduced.  A  steady-state  heat  transfer  analysis  is  also 
performed  for  the  FGM  coating/substrate  system;  the  increase  of  the  coating  thickness  owing  to 
the  functional  gradation  is  predicted. 

In  functionally  graded  coatings  made  of  layered  or  particle/matrix  type  metal-ceramic 
composites,  ligaments  of  the  metal  phase  often  bridge  the  crack  in  its  wake,  exerting  a  closure 
force  to  the  crack  surfaces,  thus  reducing  the  crack  tip  stress  intensity.  To  quantify  this  effect, 
a  crack  bridging  analysis  is  carried  out  to  predict  crack  propagation  in  coatings  made  of 
functionally  graded  materials  (FGM).  The  FGM  coating  is  taken  to  be  a  ceramic/metal 
composite  with  its  gradation  characterized  by  the  local  volume  fractions  of  metal  and  ceramic 
phases.  Fracture  in  the  FGM  coating  is  resisted  by  the  plastic  deformation  of  metal  ligaments 
in  the  crack  wake  that  bridge  the  crack;  the  crack  bridging,  however,  is  not  uniform.  A 
position-dependent  crack  bridging  model  is  developed  taking  into  account  the  coating  gradation 
and  metal  plasticity.  The  model  is  subsequently  used  in  a  finite  element  analysis  to  predict  the 
reduced  fracture  driving  force.  It  is  found  that  crack  bridging  in  the  FGM  coating  can 
significantly  reduce  the  crack  tip  stress  intensity.  It  is  also  found  that  coating  gradation  has  a 
strong  influence  on  the  fracture  driving  force  and  the  crack  length  at  arrest.  The  present  finite 
element  model  can  be  extended  readily  to  study  the  effect  of  large-scale  plastic  deformation  on 
crack  growth  in  a  FGM  coating. 

Due  to  the  performance  requirements  of  high  temperature  applications  of  functionally 
graded  coating/metal  substrate  systems,  most  of  the  FGM  coatings  are  ceramic-rich  near  the 
coating  surface.  It  is  therefore  necessary  to  analyze  fracture  of  edge  cracks  in  FGM  coatings. 
A  special  feature  pertaining  to  the  fracture  behavior  of  such  coatings  is  the  position-dependent 
crack  bridging  owing  to  the  metal  ligaments  in  the  crack  wake.  To  guide  the  design  of  coating 
gradation  so  that  cracking  damage  can  be  minimized,  a  fracture  mechanics  study  is  made  of  the 
crack  driving  force  for  bridged  edge  cracks  in  functionally  graded  ceramic/metal  coatings  using 


a  singular  integral  equations  approach.  The  metal  substrate  and  the  ceramic/metal  coating  are 
taken  as  linear  elastic,  with  the  elastic  properties  of  the  coating  varying  through  the  film 
thickness.  The  crack  density  is  assumed  to  be  relatively  small  such  that  the  crack-crack 
interactions  can  be  neglected.  Calculations  are  made  for  the  stress  intensity  factors  at  the  tip  of 
the  cracks  in  the  coating  as  determined  by  the  coating  gradation,  crack  length,  the  bridging 
characteristics,  and  the  applied  load.  It  is  found  that  crack  bridging  can  significantly  reduce  the 
crack  driving  force.  It  is  also  found  that  the  gradation  of  the  coating  affects  the  crack  driving 
force  through  an  interplay  between  the  elastic  properties  of  the  coating  and  the  crack  bridging 
characteristics.  Such  an  interplay  is  systematically  quantified. 

To  guide  the  design  of  thermal  barrier  coatings  (TBCs)  for  advanced  turbine  engines,  a 
micromechanics  analysis  is  carried  out  for  delamination  cracking  in  TBCs  consisting  of  a 
ceramic  top  coating,  a  functionally  graded  bond  coat,  and  a  metal  substrate.  The  FGM  bond 
coat  is  taken  to  be  a  multilayered  ceramic/metal  composite  with  its  thermomechanical  properties 
changing  along  the  thickness  direction  of  the  coating.  The  energy  release  rate  of  the 
delamination  crack  is  determined  analytically  as  a  function  of  the  coating  gradation,  the  location 
of  the  crack,  the  thickness  ratios  of  the  coating/substrate  system.  Finite  element  calculations 
have  also  been  performed  to  give  energy  release  rate  of  the  crack  when  the  crack  length  is 
relatively  small.  As  observed  experimentally,  cracks  perpendicular  to  the  interface  may  initiate 
in  the  top  ceramic  coating  and  propagate  into  the  bond  coat;  the  driving  force  of  such  cracks  are 
quantified  using  a  finite  element  methods.  It  is  found  that  it  is  quite  beneficial  to  use  a  FGM 
bond  coat.  However,  it  is  remain  to  be  seen  how  the  functionally  gradation  would  affect  the 
resistance  of  the  coating  to  oxidation. 

In  order  to  predict  the  fatigue  behavior  of  advanced  ceramics  and  ceramic  coatings  such 
as  alumina,  silicon  carbide  and  silicon  nitride,  a  model  is  developed  for  fatigue  growth  of  an 
edge  crack  in  nontransforming  ceramics  based  on  a  large-scale  crack  bridging  analysis.  The 
crack  bridging  is  assumed  to  be  owing  to  grain  boundary  sliding,  and  the  fatigue  behavior  of 
the  ceramics  is  modeled  using  a  frictional  wear  law  that  accounts  for  the  shielding  degradation. 
Under  cyclic  loading  conditions,  the  fatigue  problem  is  analyzed  by  solving  a  set  of  singular 
integral  equations.  It  is  found  that  the  competition  between  shielding  degradation  and  shielding 
accumulation  results  in  a  modified  Paris’  law.  It  is  also  found  that  the  dependence  on  the 
maximum  crack  tip  stress  intensity  can  be  very  strong,  but  that  on  the  R-ratio  is  only 
moderate.  Systematic  predictions  are  made  of  the  effects  of  initial  crack  length,  load  ratio,  and 
the  bridging  characteristics.  The  models  developed  can  be  further  extended  to  include  the  effect 
of  nonuniform  distribution  of  bridging  grains. 
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DELAMINATION  CRACKING  IN  FUNCTIONALLY  GRADED 
COATING/METAL  SUBSTRATE  SYSTEMS 

G.  BAO  and  H.  CAI 

Department  of  Mechanical  Engineering,  The  Johns  Hopkins  University,  Baltimore,  MD  21218,  U.S.A. 
(Received  8  May  1996;  accepted  18  June  1996) 

Abstract — A  micromechanics  analysis  is  carried  out  for  delamination  cracking  in  functionally  graded 
coating/substrate  systems.  The  FGM  coating  is  taken  to  be  a  ceramic/metal  composite  with  its 
thermomechanical  properties  changing  with  position  in  the  coating  along  the  thickness  direction.  Based 
on  a  linear  fracture  mechanics  analysis,  the  energy  release  rate  of  the  delamination  crack  is  determined 
as  a  function  of  the  coating  gradation,  the  location  of  the  crack,  the  thickness  ratio  of  the  coating/substrate 
system,  and  the  elastic  properties  of  the  ceramic  and  metal  phases.  Formulae  are  given  for  the  stress 
intensity  factors  and  the  mode  mixity  of  the  system.  Buckling  of  the  FGM  coating  due  to  compressive 
residual  stresses  is  analyzed  using  a  thin  plate  theory;  the  buckle-driven  delamination  is  modeled 
accordingly.  It  is  found  that  functional  gradation  of  the  coating  can  significantly  reduce  the  fracture 
driving  force  of  the  delamination  crack  in  both  edge-delamination  and  buckle-driven  delamination  cases. 
A  steady-state  heat  transfer  analysis  is  performed  for  the  FGM  coating/substrate  system;  the  increase  of 
the  coating  thickness  owing  to  the  functional  gradation  is  quantified.  Copyright  ©  1997  Acta  Metallurgica 
Inc. 


1.  INTRODUCTION 

A  functionally  graded  material  (FGM)  is  a  composite 
whose  composition  varies  from  place  to  place 
according  to  performance  requirements.  For 
example,  a  FGM  coating  may  have  pure  ceramic  at 
the  surface  and  pure  metal  at  the  coating/substrate 
interface.  Recent  development  of  FGM  has  demon¬ 
strated  that  such  materials  have  the  potential  to  enjoy 
a  wide  range  of  thermal  and  structural  applications, 
including  thermal  gradient  structures,  wear  and 
corrosion  resistant  coatings  and  metal/ceramic 
joining  [1,2]. 

Functional  gradation  opens  a  new  avenue  for 
optimizing  both  material  and  component  structures 
to  achieve  high  performance  and  material  efficiency; 
it  also  posts  many  challenging  mechanics  problems  in 
the  theoretical  and  experimental  determination  of  the 
microstructure-property  relationships,  and  in  corre¬ 
lating  design,  processing  and  performance  of  FGMs 

[3]. 

To  establish  the  fundamental  relationship  between 
material  gradation  and  thermomechanical  properties 
of  FGMs,  extensive  studies  have  been  carried  out  of 
the  effective  properties  [4,  5],  thermal  stress  distri¬ 
butions  [6-8]  and  fracture  in  functionally  graded 
materials  [9-13].  In  particular,  the  fracture  mech¬ 
anics  analysis  of  Erdogan  and  associates  has  been 
extended  by  Noda  and  Jin  to  include  thermal  load 
[14-16],  and  by  Bao  and  Wang  [17]  to  include  the 
effects  of  different  material  gradation.  Experimental 
observations  of  cracking  in  FGMs  have  been  made 
by  several  investigators,  including  surface  crack 


initiation  in  thermal  barrier  coatings  [18]  and 
multiple  cracking  in  a  NiAl-Al203  FGM  layer  under 
bending  [19].  The  present  work  focuses  on  the 
delamination  cracking  in  FGM  coating/substrate 
systems. 

Ceramic  coatings  have  been  developed  for  use  in 
thermal  gradient  structures,  automotive  engines  and 
cutting  and  grinding  tools  to  protect  the  surfaces 
from  melting,  wear,  corrosion  and  oxidation. 
However,  due  to  thermal  expansion  mismatch 
between  the  coating  and  substrate,  delamination 
cracking  and  coating  spalling  can  occur  upon 
experiencing  thermal  cycling  and  mechanical 
loading,  causing  premature  failure  of  the  component 
[20].  Similar  damage  mechanisms  exist  in  other 
compressed  films  [21],  Extensive  mechanics  studies 
on  delamination  and  spalling  in  film/substrate 
systems,  with  both  the  film  and  substrate  being 
homogeneous,  have  been  carried  out  [22-25].  A 
review  of  the  mechanics  models  for  delamination 
cracking  in  such  systems  can  be  found  in  Hutchinson 
and  Suo  [26], 

To  alleviate  the  problems  encountered  in  using  the 
pure  ceramic  coatings,  functionally  graded  materials 
have  been  developed  for  use  as  thermal  barrier 
coatings  (TBC)  and  as  protective  coatings  for 
hard-metal  tools  [27],  FGM  coatings  have  the 
potential  to  simultaneously  reduce  thermal  expansion 
mismatch,  increase  interface  bonding  strength,  and 
enhance  coating  toughness.  To  realize  this  potential, 
an  effort  is  made  in  the  present  study  to  establish  the 
relationship  between  delamination  cracking  of  the 
coating  and  coating  gradation,  aiming  to  provide 
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Fig.  1.  A  schematic  of  delamination  cracking  in  a 
functionally  graded  coating/metal-substrate  system  under 
applied  forces  and  moments. 


guidance  to  both  the  coating  design  and  the  life 
prediction  of  the  FGM  coating/substrate  systems. 

2.  PRELIMINARIES 

Consider  a  functionally  graded  ceramic/meta! 
coating  on  a  homogeneous,  isotropic  metal  substrate, 
as  shown  schematically  in  Fig.  1 .  The  coating  is  taken 
to  be  a  ceramic  of  an  inclusion/matrix  type  shown  in 
Fig.  2(a),  or  a  functionally  graded  multilayer  depicted 
in  Fig.  2(b).  Based  on  the  performance  requirements 
for  thermal  barrier  coatings  and  coatings  in 
wear-related  applications,  in  this  study,  only  coatings 
that  are  ceramic-rich  near  the  coating  surface  and 
metal-rich  near  the  coating/substrate  interface  are 
considered.  For  convenience,  the  y-axis  is  set  along 
the  coating  thickness  direction,  the  x-axis  lies  within 
the  coating/substrate  interface,  as  illustrated  in  Fig.  1. 
At  any  position  y  in  the  ceramic/metal  coating,  the 
local  volume  fraction  of  ceramic  is  assumed  to  be  g(y) 
which  can  be  used  to  characterize  the  coating 
gradation.  Generally  speaking,  g(y )  can  be  any 
non-singular,  non-negative  function  of  y.  To  gain 
insight,  it  is  assumed  that  g(y)  obeys  a  pure 
power-law  type  relation 

g(y)  =  (y/hy  0) 

where  h  is  the  coating  thickness,  and  m  is  a  material 
parameter.  The  gradation  given  in  equation  (1) 
implies  that  the  coating  always  has  100%  ceramic  at 


the  surface  (i.e.  g(h)  =  1)  and  100%  metal  at  the 
interface.  Though  idealized,  the  coating  gradation 
defined  in  equation  (1)  can  be  used  to  model 
approximately  the  FGM  coatings  of  both  the 
inclusion/matrix  type  shown  in  Fig.  2(a),  and  the 
multilayered  type  shown  in  Fig.  2(b).  The  total 
volume  fraction  /c  of  ceramic  in  the  FGM  coating  is 
related  to  the  exponent  m  in  equation  (1)  by 

fc  =  1  Km  +  1).  (2) 

The  coating  gradation  given  in  equation  (1)  is  similar 
to  that  used  by  Drake  et  al.  [7]  and  Bao  and  Wang 
[17]  in  their  studies  of  FGM  systems. 

Due  to  the  heterogeneous  nature  of  a  functionally 
graded  material,  it  is  quite  difficult  to  obtain  the  exact 
solution  of  the  effective  elastic  moduli  and  thermal 
expansion  coefficient  of  a  FGM  coating  in  terms  of 
the  relative  volume  fractions  and  thermomechanical 
properties  of  the  metal  and  ceramic  phases.  For  a 
FGM  coating  in  the  inclusion/matrix  type  (Fig.  2(a)), 
the  existing  micromechanical  cell  models  are  no 
longer  valid  since  no  unit  cell  containing  a  single 
particle  can  be  used  to  represent  the  whole  composite 
body.  For  a  FGM  with  a  layered  structure  as  that 
depicted  in  Fig.  2(b),  the  existing  micromechanical 
cell  models  are  valid  only  if  the  thickness  of  each 
individual  layer  is  much  larger  than  the  average 
particle  size. 

As  a  first-order  approximation,  in  this  study 
Young’s  modulus  E{y)  and  thermal  expansion 
coefficient  a(y)  of  the  FGM  coating  are  assumed  to 
have  power-law  type  expressions 

E(y)  =  Em  +  {E^,Em)(^j  (3a) 

a  O')  =  am  +  (ac  -  am)^J  (3b) 

where  Em  and  £c  are  Young’s  moduli,  and  am  and  ac 
the  thermal  expansion  coefficients  of  the  metal  and 
ceramic  phases,  respectively,  and  n  and  r  are 
exponents  determined  by  the  coating  gradation. 


Ceramic 

Matrix 


Particles 


Ceramic-rich 

Layer 


Metal-rich 

Layer 


(a) 


(b) 


Fig.  2.  Typical  microstructures  of  functionally  graded  coatings;  (a)  an  inclusion/matrix  structure;  (b)  a 

multilayered  structure. 
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Fig.  3.  Young’s  modulus  of  the  coating  as  a  function  of 
position  as  determined  using  a  power-law  type  model  with 
different  gradation  exponent. 


Shown  in  Fig.  3  are  curves  for  E(y)/Em  as  a  function 
of  y/h  for  Ec/Em  =  3  for  n  =  0.5,  1  and  2.  Note  that 
n  —  r  =  0  corresponds  to  a  pure  ceramic  coating.  For 
simplicity,  in  the  case  studies,  it  is  assumed  that 
Poisson’s  ratios  of  the  metal  and  ceramic  are  the 
same,  i.e. 


Vm  =  v„  =  V.  (4) 

Practically,  the  exponents  n  and  r  in  equations  (3a) 
and  (3b)  can  be  determined  from  measurements  of 
the  position-dependent  Young’s  modulus  and  the 
thermal  expansion  coefficient  of  the  FGM  coating. 
They  can  also  be  obtained  approximately  from  the 
predictions  of  Young’s  modulus  and  the  thermal 
expansion  coefficient  of  the  FGM  coating  using  a 
micromechanics  cell  model  [28,  29]  in  conjunction 
with  the  coating  gradation  (1).  Although  Poisson’s 
ratios  of  the  metal  and  ceramic  phases  are  usually 
different,  and  Poisson’s  ratio  of  the  FGM  coating 
varies  with  position,  the  effect  of  using  equation  (4) 
as  an  approximation  is  expected  to  be  small. 


3.  EDGE  DELAMINATION 

As  depicted  in  Fig.  1,  the  FGM  coating/metal 
substrate  system  containing  an  edge  delamination 
crack  is  assumed  to  be  in  plane  strain  conditions  and 
subjected  to  axial  loads  (force  per  unit  width) 
P,,P2,  P3  and  bending  moments  (moment  per  unit 
width)  Mi,  M2  and  M3.  Both  the  FGM  coating  and 
the  metal  substrate  are  taken  to  be  isotropic,  linear 
elastic,  with  Young’s  moduli  E(y),  Em  and  Poisson’s 
ratio  v(y),  vm  for  the  coating  and  substrate,  respect¬ 
ively;  plasticity  of  the  metal  phase  and  the  possible 
crack  bridging  mechanisms  are  neglected.  The  crack, 
which  is  located  a  distance  hi  from  the  coating 
surface,  has  length  a  much  larger  than  the  coating 
thickness  h.  Attention  in  this  study  is  restricted  to 
situations  where  hi  <  h,  although  the  same  approach 
can  be  used  for  hi  >  h.  When  hi  =  h,  the  delamina¬ 
tion  crack  lies  within  the  interface  between  the 


coating  and  the  substrate.  Note  that  the  thermomech¬ 
anical  properties  of  the  coating  defined  in  equations 
(3)  and  (4)  have  no  discontinuities  within  the  coating 
nor  at  the  interface  y  =  0.  The  delamination  crack, 
therefore,  can  be  analyzed  without  introducing  the 
complex  stress  intensity  factor  [30],  However,  the 
fracture  toughness  of  the  interface  may  still  be  lower 
than  that  of  the  metal  substrate  owing  to  the 
deposition  process. 

3.1.  General  results 

The  energy  release  ^  of  the  delamination  crack  can 
be  derived  using  an  elementary  beam  theory  [25]. 
Assuming  that  the  three  edges  have  linear  strain 
distributions,  the  positions  of  the  neutral  axes  5i,  d2 
and  5  (Fig.  1)  can  be  determined  by 

<5,  =  f  -  (h  -  hi)  (5) 

/-O 

c  H  2(o>i  —  AQ  +  H(wo  —  /q)  , 

2  "  2  +  2 (EmH  +  (o0-  Ac)  W 

c  _  H  2(0 1  +  H(Oq  .r,s 

2  +  2(EmH  +  (o0)  K) 

where  H  is  the  thickness  of  the  substrate,  Em  =  Emj 
(1  —  v2)  is  the  plane  strain  Young’s  modulus  of  the 
metal  substrate, 


n 


£00  dy 

7i  -  h  i 

(8a) 

n 

E(y)y  d y 

It -In 

(8b) 

|  Eiy)dy 

(9a) 

Jo 


0)1  = 


E(y)y  d y 
Jo 


(9b) 


where  E(y )  =  E(y)/[  1  —  v2(y)]  is  the  plane  strain 
Young’s  modulus  of  the  FGM  coating.  Since  a»h, 
the  energy  release  rate  is  simply  the  difference  of  the 
strain  energy  per  unit  length  per  unit  width  stored  in 
the  edges  far  behind  and  far  ahead  of  the  crack  tip 


.  _  P\  .  M\ 

220  +  2(A2  -  Af/Ao) 


+ 


2(co0  —  Ao  *F  HE m) 


Aft 


Ml 


2(4>2  +  (j>  m)  2(a>o  +  HEm)  2{\[ii  +  if/m) 


(10) 


where 


A2  — 


'h 

E{y)y 2  d y 

J/j  -  h\ 


(11  a) 


(jh  =  (//  -  <52)2(co„  -  A„)  +  2 (H  -  52){a>i  -  ki) 

+  CO  2  —  ^2  (lit)) 
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i fa  =  (H-  Sfco o  +  2 (H  -  6)0),  +  a>2  (1  lc) 

with 


%h 

w2  =  E(j>)y 2  dy 

Jo 

(lid) 

and 

^  =  ~[(ff  -  s2y  +  si] 

(12a) 

iltm  =  ^[(H-sy  +  n 

(12b) 

The  expression  of  ^  given  in  equation  (10)  can  be 
further  simplified  using  a  superposition  scheme  [25], 

„  1  [  P2  M2  2  PM  . 

rn  +  wn  +  . _  sm  y 

2Em  \_Uh  Vh  /[JVh2 

(13) 

where 

P  =  P,  -  C,P3  -  C2M3/h 

(14a) 

M  =  M,  -  M3C3 

(14b) 

with 

1  Wo  +  HEm 

(15a) 

h[X,  +  (H-  <5)2»] 

02  +  0m 

(15b) 

2.2  —  X2/Ao 

3  02  +  0m 

(15c) 

The  non-dimensional  parameters  U,  V  and  sin  y  in 
(13)  are  given  by 

I-P 

, fl  1  .  (H  —  S2  +  2i/ilo)2"l 

U  m  |_20  w0  -  20  +  HEm  +  <p2  +  0m  J 

(16a) 


of  the  FGM  coating  and  the  substrate.  Once  w  is 
given,  the  stress  intensity  factors,  and  the  phase  angle 

0 

0  =  tan-1^^  (18) 

can  be  fully  determined. 

3.2.  Total  energy  release  rate 

The  results  given  in  equations  (5)— ( 1 7)  are  valid  for 
a  FGM  coating  with  any  elastic  moduli  E(y)  and 
v(y).  To  gain  insight,  in  the  following,  we  assume  that 
E(y )  has  a  power-law  form,  equation  (3a), 

E{y)  =  Em  +  (E,  -  Em)(^j  (3a) 

and  that  Poisson’s  ratio  of  the  FGM  coating  is 
constant  (  =  vm).  Consequently,  the  parameters 
2o,  Wo,  (0,  and  co2  defined  in  equations  (8),  (9) 
and  (11)  can  be  expressed  in  a  closed-form;  they  are 
for  i  =  1,  2  and  3, 


The  strain  energy  release  rate  ^  of  the  delamina¬ 
tion  crack  is  a  function  of  the  applied  loads,  and  the 
geometric  and  material  parameters.  In  a  non-dimen¬ 
sional  form,  we  have 


—  =  F  IP 

y  -  k^m'1 


$E m  ,/h,  HEC 


sin  y  Emh2(H  —  62  +  2,/20)  n  ^ 
<h  +  4> m  '  (  ; 

It  is  readily  shown  that  for  a  homogeneous  coating, 
the  expressions  in  equation  (10),  (15)  and  (16)  are 
identical  to  those  given  in  Suo  [25], 

The  delamination  crack  shown  in  Fig.  1  is,  of 
course,  a  mixed-mode  crack.  The  stress  intensity 
factors  K,  and  Kn  can  be  derived  from  equation  (13) 
based  on  linearity  and  dimensionality. 


Kt  -  —  cos  co  -1 - - - sin(o>  +  y)  (17a) 

Jtmj  JlWV 

Kn  =  — — —  sin  co  —  — ^ . cos(w  +  y)  (17b) 

j2hU  jWV 

where  co  is  a  function  of  the  thickness  ratio  hjH,  the 
location  of  the  crack  hjh,  and  the  material  properties 


where  0  is  a  non-dimensional  function,  a  is  a  stress 
measure  related  to  the  applied  loads  but  not 
dependent  on  the  parameters  in  0.  For  any  FGM 
coating  on  a  homogeneous  substrate  subjected  to 
applied  loads  P,,P2,  P2  and  M,,M2,  M3  as  shown  in 


Fig.  4.  A  schematic  of  delamination  cracking  in  a 
functionally  graded  coating/metal-substrate  system  under 
far-field  uniform  straining. 
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Fig.  5.  (a)  The  normalized  energy  release  rates  3  as  a  function  of  the  thickness  ratio  A/Lf  for  a  crack  with 
location  Ai/A  =  0.4,  0.6,  0.8  and  1.0;  (b)  ^  as  a  function  of  the  crack  location  ht/h  for  the  gradation 
exponent  n  =  0.5, 1.0  and  2.0  for  h)H  =  0.01. 


i 


Fig.  1,  (j>  can  be  calculated  readily  using  the  above 
equations. 

To  illustrate  the  trends  in  the  normalized  &  in 
equation  (20),  consider  the  FGM  coating/metal 
substrate  system  subjected  to  a  uniform  applied  strain 
£0,  as  shown  in  Fig.  4.  We  have  in  this  case 

P  =  £o2o,  M  =  0  (21) 

and  the  stress  a  in  (20)  can  be  taken  as  a  =  eaEm. 
Displayed  in  Fig.  5(a)  is  the  normalized  energy 


release  rate  as  a  function  of  the  thickness  ratio  hjH 
for  EJEm  —  3,  n  =  1  for  hi/h  =  0.4,  0.6,  0.8  and  1.0. 
It  is  seen  that  the  total  strain  energy  release  rate  ^ 
increases  with  increasing  thickness  ratio  h/H. 
Further,  with  a  fixed  h/H,  <3  increases  with  increasing 
hi/h.  Clearly,  the  location  of  the  delamination  crack 
depends  on  3  as  well  as  the  toughness  E  of  the  FGM 
layer  which  is  also  a  function  of  position.  When  the 
FGM  coating  is  ceramic-rich  at  the  surface  and 
metal-rich  at  the  coating/substrate  interface,  F  is 
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expected  to  be  higher  near  the  interface.  Conse¬ 
quently,  the  delamination  crack  is  likely  to  be  located 
some  distance  away  from  the  interface.  To  show  the 
effect  of  material  gradation,  in  Fig.  5(b),  the 
normalized  ^  is  plotted  against  the  relative  location 
of  the  crack  hfh  for  Ec/Em  =  ?i,  h/H  =  0.01  for 
n  =  0.5,  1  and  2.  The  fracture  driving  force  is  lower 
when  n  is  larger,  as  can  be  seen  from  Fig.  5(b). 

3.3.  Stress  intensity  factors 

The  stress  intensity  factors  Ki  and  Kn  of  the 
delamination  crack  can  be  calculated  using  equation 
(17)  if  the  parameter  co  is  known.  When  the  gradation 
of  the  FGM  coating  obeys  equation  (3a),  co  is  a 
function  of  the  thickness  ratio  H/h,  crack  location 
hfh,  and  material  parameters  Ec/Em  and  n,  i.e. 


H  E 
h  '  Em 


(22) 


To  quantify  co,  a  finite  element  analysis  was  carried 
out  using  the  commercial  code  ABAQUS  for  the 
FGM  coating/metal  substrate  system  depicted  in 
Fig.  1 .  The  FGM  coating  is  divided  into  200  layers, 
the  first  layer  in  the  coating  away  from  the  interface 
is  assumed  to  have  the  same  Young’s  modulus  as  of 
the  metal  matrix.  For  «  =  1,2,...,  99,  the  (2«)th 
and  (2 n  +  l)th  layers  have  the  same  Young’s 
modulus  E(y„)  with  y„  =  0.01  nh.  The  delamination 
crack  in  the  FGM  coating  is  assumed  to  be  located 
in-between  (2n)th  and  (2 n  +  l)th  layers;  compli¬ 
cations  due  to  the  oscillatory  singular  crack-tip  fields 
of  an  interface  crack  is  thus  avoided.  Shown  in  Fig.  6 
is  an  example  of  the  finite  element  mesh  for  the  FGM 
coating.  Clearly,  in  this  finite  element  model  the 
gradation  is  not  continuous,  nor  is  it  in  the  real 
material.  This  model  mimics  the  microstructure  of  a 
fine  multilayered  FGM  coating,  or  an  inclusion/ 
matrix  type  FGM  coating  with  gradation  E(y)  given 
in  equation  (3a).  The  stress  intensity  factors  Kt  and 
Ku  were  determined  through  the  calculated  crack 
opening  displacements  [33].  The  accuracy  of  the 
calculated  K,  and  K„  was  assessed  by  taking  the 
coating  as  homogeneous  and  comparing  the  resulting 
/Gvalues  with  that  obtained  by  Suo  [25]. 


FGM  Coating 


Crack-tip 

Region 

Fig.  6.  A  typical  finite  element  mesh  used  in  calculating  the 
stress  intensity  factors  of  a  delamination  crack  in  the 
functionally  graded  coating. 


Crack  Location,  h,/h 


coating/metal  substrate  system  with  n  =  1,  h/H  =  0.1. 


Shown  in  Fig.  7(a)  is  co  as  a  function  of  the  relative 
location  of  the  crack  ht/h  for  H/h  =  10,  Ec/Em  =  3, 
n  =  1 .  With  less  than  1  %  error,  the  curve  for  co  shown 
in  Fig.  7(a)  can  be  represented  by 

co  =  51.308  +  6.688(/z,//t) 

+  13.038(/i,//t)2  -  16.301(/t,//03.  (23) 

The  corresponding  phase  angle  of  the  delamination 
crack,  i/r  =  tan~'(Kn/Ki)  is  plotted  in  Fig.  7(b).  It  is 
interesting  to  note  that  the  minimum  value  of  t/r 
occurs  at  hfh  =  0.8  which  implies  that,  if  the  coating 
has  a  constant  toughness  and  if  fracture  is  locally 
mode-I  controlled,  the  delamination  crack  will  be  at 
hi  =  0.8/r . 

As  a  limiting  case,  the  delamination  crack  located 
at  the  coating/substrate  interface  (i.e.  hi  =  h)  is 
analyzed  using  a  finite  element  method  and  the 
corresponding  function  co  is  obtained.  Displayed  in 
Fig.  8(a)  is  co  as  a  function  of  the  thickness  ratio 
rj  =  hjH  for  n  =  0.5,  1.0  and  2.0.  The  points  are  the 
finite  element  results;  the  solid  lines  are  the 
curve-fittings.  It  is  seen  that  for  different  values  of 
n,  co  can  be  represented  approximately  as  a  linear 
function  of  hjH.  Shown  in  Fig.  8(b)  are  co  vs  \/n 
curves  for  Hjh  =  1.0,  2.0,  3,0,  4.0,  5.0  and  10.0. 
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Within  a  few  percent  error,  these  curves  can  be 
reproduced  by 

co  =  coo («)  -  (3  +  1  ln)(h/H)  (24a) 
co0(n)  =  52.1  -  7.3 1/«  +  15.264/«2  -  4.894/n3. 

(24b) 

The  stress  intensity  factors  for  a  crack  located  at  the 
coating/substrate  interface  can  be  calculated  using 
equations  (17)  and  (24). 

4.  BUCKLE-DRIVEN  DELAMINATION 

Delamination  cracking  often  occurs  at  the  coating/ 
substrate  interface  when  the  coating  buckles  due  to 
the  compressive  thermal  residual  stress  [26].  This 
buckle-driven  delamination  cracking  can  also  occur 
in  a  FGM  coating/metal  substrate  system,  although 
the  residual  stress  is  somewhat  reduced.  A  new 
feature  for  the  FGM  system  is  that  the  delamination 
crack  may  be  located  in  the  coating  away  from  the 
interface,  since  at  the  interface  the  residual  stress  is 
lower  and  the  toughness  can  be  much  higher  than 


Thickness  Ratio,  T|  =  h/H 


Fig.  8.  (a)  The  parameter  co  as  a  function  of  the  thickness 
ratio  h/H  for  the  gradation  exponent  n  =  0.5,  1.0  and  2.0; 
(b)  the  parameter  a  as  a  function  of  the  coating  gradation 
exponent  n  for  thickness  ratio  Hjh  =  1.0,  2.0,  3.0,  4.0,  5.0, 
and  10.0.  The  delamination  crack  is  located  at  the 
coating/substrate  interface. 


Local  Loading 
of  Crack 


Buckled 


Fig.  9.  A  schematic  of  a  buckle-driven  delamination  crack 
in  a  FGM  coating  on  an  infinitely  deep  metal  substrate. 


that  at  other  places  in  the  FGM  coating.  Attention 
in  this  paper  is  restricted  to  the  case  in  which  the 
thickness  of  the  substrate  is  much  larger  than  that  of 
the  coating,  i.e.  h/H^>0.  This  will  greatly  simplify  the 
buckling  analysis  of  the  film. 

Consider  an  one-dimensional  blister  shown  in 
Fig.  9.  The  FGM  coating  is  assumed  to  contain  a 
debonded  region  of  width  2b  located  a  distance  h\ 
from  the  coating  surface.  Under  an  equi-biaxial 
compressive  in-plane  stress  axx  =  ov  =  a  (e.g.  a 
thermal  residual  stress),  the  debonded  region  may 
buckle  up,  inducing  a  plane  strain,  mixed  mode 
delamination  crack.  The  energy  release  rate  ^  of  the 
delamination  crack  can  be  obtained  by  analyzing  the 
local  loads  P  and  M  on  the  buckled  film,  then 
inserting  them  into  equation  (13). 

To  calculate  the  loading  of  the  delamination  crack, 
it  is  assumed  that  h\«b  so  that  the  debonded  region 
shown  in  Fig.  9  can  be  represented  by  a  wide, 
clamped  Eular  column  of  width  2b.  Following 
Hutchinson  and  Suo  [26],  the  displacement  v  of  the 
film  in  the  y-direction  is  governed  by 


=r  d4t>  ...  d2!)  A 

D  dx4  {p  N r)  dx2  ~  ° 


(25) 


where  P  is  the  change  in  the  resultant  force  in  film  due 
to  buckling,  N,  is  the  resultant  force  owing  to  the 
thermal  residual  stress,  and  D  the  bending  stiffness 


D  = 


(y  —  h  +  h,  —  d,)2E(y)  dy  =  h  -  k\jk0. 


(26) 

Solving  equation  (25)  with  the  boundary  conditions 
del 


v\x  =  ±b  =  0, 


dx 


=  0 


(27) 


and  taking  the  lowest  non-zero  eigen-value,  we  have 


v{x)  =  \m 


i  ,  nx 

1  4-  cos  ~r 


(28a) 
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where  [26] 


£ 


K0) 

h 


4[D(N, 
h  \NC 


(28b) 


with 


(28c) 


The  loading  P  and  M  can  be  given  readily  from  the 
above  solution 


P  =  N,  -  Nc  (29a) 


M  = 


n2Dh  K 


(29b) 


The  energy  release  rate  of  the  buckle-driven 
delamination  crack  is  then  obtained  from  (13) 


Temperature  Drop,  AT  /ATc(=Nr/Nc) 


Fig.  10.  The  normalized  energy  release  rate  of  a 
buckle-driven  delamination  crack  as  a  function  of  the 
temperature  drop  for  crack  locations  h\/h  =  0.2,  0.4,  0.6,  0.8 
and  1 .0.  Note  that  the  energy  release  rates  for  hi/h  —  0.4  and 
lh/h  =  1.0  are  almost  identical. 


K 

220 


temperature  drop  AT/ATQ  =  Nr/Nc  for  Ec/Em  =  3, 
(30)  r  =  n  =  1  for  E/h  =  0.2,  0.4,  0.6,  0.8  and  1.0.  The 
reference  energy  rate  is  given  by 


The  thermal  residual  stress  distribution  in  a  FGM 
system  can  be  calculated  readily  once  the  thermo¬ 
mechanical  properties  of  the  coating  and  the 
substrate  are  given  [31].  As  mentioned  earlier,  in  this 
study  we  assume  that  the  variation  of  both  Young’s 
modulus  E(y)  and  the  thermal  expansion  coefficient 
a(y)  in  the  FGM  coating  obeys  a  power-law  type 
relation  defined  by  equation  (3).  Further,  Poisson’s 
ratio  of  the  coating  is  taken  to  be  constant  (  =  vm). 
Based  on  these  assumptions,  the  residual  stress 
distribution  in  the  FGM  coating  is  given  by 

ffroo  =  £(y)W««  ~aM  pi) 

1  vm 


Em(am  -  QA  T 
1  -  vm 


(33) 


The  critical  temperature  drop  A  7)  is  the  temperature 
drop  at  which  the  film  buckles 


TvP{\  -  Vm) 

c  Fra(am  -  ac)hb2 


+ 


_ 

n  +  r  +  1 


where  AT  is  the  temperature  change.  Typically 
AT<0  and,  for  a  FGM  metal/ceramic  coating, 
a(y)  am;  a t(y)  in  equation  (31)  is  therefore 
compressive.  The  resultant  force  N,  owing  to  the 
thermal  residual  stress  is 


Nt  = 


'it 

ov( v)  d y 

4  h  —  h  | 


Fm(am  —  ac)hAT 
1  -  vm 


1 

r+  1 


1  f  Ec 
+  n  +  r  +  1  \Em 


(32) 

Plotted  in  Fig.  10  are  curves  of  the  normalized 
energy  release  rate  against  non-dimensional 


It  is  noteworthy  that  when  a  crack  is  located  at 
h\jh  =  0.6  the  fracture  driving  force  is  higher  than 


Fig.  1 1 .  The  nondimensional  critical  temperature  drop  as  a 
function  of  the  gradation  exponent  n  for  a  buckle-driven 
delamination  crack  with  location  hi/h  =  0.4,  0.6,  0.8  and 
1.0. 
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Temperature  Drop,AT/ATc(=N,/Nc) 

Fig.  12.  The  normalized  energy  release  rate  of  a 
buckle-driven  delamination  crack  as  a  function  of  the 
temperature  drop  for  coating  gradation  exponent 
n  =  0, 0.5,  1.0  and  2.0.  The  crack  is  located  at  the 
coating/substrate  interface.  Note  that  n  =  0  corresponds  to 
a  pure  ceramic  coating. 


that  with  other  crack  locations.  The  curve  with 
/ti//j  =  0.4  is  almost  identical  to  that  with  h,/h  =  1.0. 
For  each  curve,  the  peak  value  of  occurs  at 
AT/ATC  =  3.0. 

The  critical  temperature  drop  A  Tc  is  an  important 
design  parameter  for  if  AT  <  A Tc,  no  buckling  of  the 
film  will  occur.  Shown  in  Fig.  1 1  is  the  non-dimen¬ 
sional  temperature  drop  as  a  function  of  the 
gradation  exponent  n  for  Ec/En  =  3,  n  =  r  =  1  for 
h]/h  =  0.4,  0.6,  0.8  and  1.0.  It  is  seen  that  the  change 
of  ATC  with  n  is  almost  linear.  Since  increasing  n 
implies  to  increase  the  total  volume  fraction  of  metal 
in  the  FGM  coating,  the  trend  in  ATC  vs  n  curves 
shown  in  Fig.  1 1  might  be  expected.  Note  also  that 
when  n  =  0,  the  FGM  coating  becomes  a  pure 
ceramic  coating,  so  the  critical  temperature  drop  is 
the  smallest  for  all  values  of  h^h.  When  the 
delamination  crack  is  located  near  the  interface,  the 
film  is  thicker  than  otherwise;  A  Tc  is  thus  higher.  The 
opposite  is  true  when  the  crack  is  located  near  the 
coating  surface.  Clearly,  for  a  fixed  temperature  drop 
(e.g.  from  the  processing  temperature  to  room 
temperature),  initial  flaws  in  the  FGM  coating 
located  near  the  coating  surface  are  more  likely  to 
cause  local  buckling. 

The  use  of  a  FGM  ceramic/metal  coating  can 
significantly  decrease  the  residual  compressive 
stresses  in  the  coating,  thus  reducing  the  possibility  of 
buckle-driven  delamination.  To  illustrate,  in  Fig.  12, 
the  normalized  energy  release  rate  is  plotted 
against  non-dimensional  temperature  drop  AT/ATQ 
(=N,/NC)  for  EJEm  =  3,  r  =  n,  h\  =  h  for 
n  =  0.0,  0.5,  1.0  and  2.0.  Note  that  with  increasing  n, 
the  energy  release  rate  decreases,  since  the  volume 
fraction  of  metal  in  the  coating  increases.  For  a  pure 
ceramic  coating  ( n  =  0),  the  fracture  driving  force  is 
much  higher  than  that  for  a  FGM  coating,  as  can  be 
seen  from  Fig.  12. 


5.  STEADY-STATE  HEAT  TRANSFER 

As  shown  clearly  in  the  previous  section,  with  a 
FGM  coating  the  buckle-driven  delamination  crack 
can  have  a  much  lower  fracture  driving  force 
compared  with  a  pure  ceramic  coating.  Flowever,  the 
FGM  coating  contains  metal  phase  which  may  have 
a  much  higher  thermal  conductivity.  Consequently, 
when  the  FGM  coating  is  used  as  a  thermal  barrier 
coating,  to  provide  the  same  degree  of  protection  to 
the  metal  substrate,  a  thicker  coating  has  to  be  used. 

To  predict  the  thickness  of  a  FGM  coating  as 
determined  by  the  heat  shielding  requirements  and 
the  coating  gradation,  consider  the  FGM  coating/ 
metal  substrate  system  shown  in  Fig.  13.  It  is  assumed 
that  the  working  temperature  (e.g.  the  temperature  of 
gas  in  the  engine)  on  the  coating  surface  is  To,  and  the 
temperature  the  metal  substrate  can  sustain  is  Tm.  For 
simplicity,  T0  is  assumed  to  be  constant;  for  a 
conservative  design  it  can  be  taken  as  the  maximum 
gas  temperature  at  coating  surface.  Further,  the 
coating  thickness  h  is  taken  to  be  much  smaller  than 
other  dimensions  of  the  system  so  that  the  edge 
effects  can  be  neglected.  The  resulting  one-dimen¬ 
sional  steady-state  heat  transfer  problem  is  governed 
by 


&=~k{y)%  (35) 

where  T(y)  is  the  temperature  distribution  in  the 
coating,  Q  is  the  heat  flux  per  unit  area,  and  kiy)  is 
the  thermal  conductivity  of  the  FGM  coating. 
Integration  of  equation  (35)  gives 


d y 

lo 


(36) 


For  a  pure  ceramic  coating  with  the  same  surface 
temperatures  and  heat  flux,  we  have 

T0  —  Tm=  (37) 

where  hc  is  the  thickness  and  kc  the  thermal 
conductivity  of  the  ceramic  coating.  Combining 


,  y 

To 

FGM  Coating 

10 

Tm 

Substrate 

Fig.  13.  A  schematic  of  the  steady-state  heat  transfer 
problem  for  a  FGM  coating/metal  substrate  system. 
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equation  (36)  and  (37)  leads  to 


'*  d s 
k(s) 


(38) 


where  s=yjh.  It  is  interesting  to  note  that  the 
thickness  ratio  /? //?c  only  depends  on  the  thermal 
conductivities  of  the  FGM  coating  and  the  ceramic. 
Obviously,  for  a  pure  ceramic  coating,  k(s)  =  kc, 
h  =  /;c;  for  a  ceramic/metal  FGM  coating,  k(s)  >  kc, 
h  >  hc  which  implies  that  when  a  thermal  barrier 
coating  is  made  of  functionally  graded  material,  a 
thicker  coating  needs  to  be  used. 

The  thermal  conductivity  k  of  a  dual-phase 
composite  with  spherical  inclusions  uniformly 
distributed  in  the  matrix  is  given  by  [32] 


k  =  km 


1  + 


(1  -  c)/ 3  +  km/(ki  -  km) 


(39) 


where  the  subscripts  i  and  m  represent  the  inclusion 
and  matrix  phases,  respectively,  and  c  is  the  volume 
fraction  of  the  inclusion.  For  a  FGM  coating, 
however,  the  definitions  of  “matrix”  and  “inclusion” 
are  position  dependent:  near  the  surface  of  the 
coating  the  ceramic  phase  is  the  matrix  and  the  metal 
phase  is  in  the  form  of  inclusion,  while  near  the 
coating/substrate  interface,  the  opposite  is  true.  At 
the  middle  of  the  coating  (the  percolation  area)  the 
ceramic  and  metal  phases  are  both  continuous; 
matrix  and  inclusion  cannot  be  well  defined.  Further, 
for  a  FGM  coating  the  inclusion  phase,  be  it  ceramic 
or  metal,  is  not  uniformly  distributed  in  the  matrix. 
The  expression  given  in  equation  (39)  is  a  good 
approximation  only  for  a  multi-layered  coating  with 
inclusion  size  much  smaller  than  the  thickness  of  each 
layer. 

To  gain  insight,  in  this  study,  we  assume  that  up 
to  the  percolation  limit  equation  (39)  can  still  be  used 
to  obtain  the  effective  thermal  conductivity  k  of  the 
FGM  coating.  The  local  volume  fraction  g  of  the 
ceramic  phase,  however,  is  a  function  of  the  position 
y  in  the  coating.  We  further  assume  that  at  the 
percolation  limit  gP,  there  is  a  sudden  change  in  the 
effective  thermal  conductivity,  as  shown  in  Fig.  14. 
Specifically,  when  g(y )  <  gp,  the  metal  phase  is  the 
matrix, 


k(y)  =  ta, 


1  + 


g  00 


[1  +  g(y)]/3  +  ta/(ta  —  km) 


■  (40) 


When  g(y)  >  gp,  the  ceramic  phase  becomes  the 
matrix, 


Fig.  14.  The  effective  thermal  conductivity  of  the  FGM 
coating  as  a  function  of  the  position  in  the  coating  for 
km/ke  =  10.  The  local  volume  fraction  of  ceramic  in  the 
coating  is  characterized  by  g (y)  =  (y/h)m.  The  solid  lines  are 
the  predictions  obtained  using  the  composite  spheres  model; 
the  dashed  line  are  upper  and  lower  bounds  based  on  the 
rule  of  mixtures.  Three  percolation  limits  (gp  =  0.4,  0.5  and 
0.6)  are  considered. 


phases.  To  simplify  the  analysis,  three  values  of  gp  are 
considered:  gp  =  0.4,  0.5  and  0.6.  Displayed  in  Fig.  14 
is  the  thermal  conductivity  of  the  FGM  coating  as  a 
function  of  the  position  y  in  the  coating  correspond¬ 
ing  to  a  pure  power-law  gradation 

g(y)  =  iy/hT.  (1) 

Also  shown  in  Fig.  14  are  the  upper  and  lower 
bounds  of  k(y)  obtained  using  the  rule  of  mixtures, 

tapper  =  gOOta  +  [1  -  g{y)]km 

fclower  =  {g(y)/kc  +  [1  -  gOOl/ta,}'1.  (42) 

Ftaving  obtained  the  effective  thermal  conductivity 
of  the  FGM  coating,  the  thickness  ratio  h/hc  can  be 
predicted  readily  using  equations  (1)  and  (38)  in 
conjunction  with  equation  (40)  and  (41).  Plotted  in 


k(y)  =  ta 


1  -g(y) 

[g<j’)V 3  +  ta/(ta  -  ta) 


(41) 


The  percolation  limit  gp  of  a  metal/ceramic 
composite  is  typically  ranging  from  0.4  to  0.6, 
depending  on  the  specific  microgeometry  of  the 


Gradation  Exponent,  m 

Fig.  15.  Predictions  of  the  coating  thickness  ratio  h/hc  as  a 
function  of  the  coating  gradation  exponent  m  based  on  the 
effective  thermal  conductivity  of  the  FGM  coating  obtained 
using  different  approaches.  The  gp  =  0.5  curve  provides  a 
good  estimate  of  the  coating  thickness  as  compared  with 
that  of  a  pure  ceramic  coating. 
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Fig.  15  is  the  thickness  ratio  h/hc  as  a  function  of  the 
coating  gradation  exponent  m  for  km/kc  =10  for 
gp  =  0.4,  0.5  and  0.6.  Shown  also  are  h/lu  vs  m  curves 
corresponding  to  the  upper  and  lower  bounds  of  the 
effective  thermal  conductivity  of  the  coating  given  in 
equation  (42).  Indeed  the  metal  phase  in  the  FGM 
coating  could  result  in  a  significantly  thicker  coating 
compared  with  the  pure  ceramic  coating.  Note  that 
the  coating  gradation  exponent  m  is  related  to  the 
total  volume  fraction  of  metal  /m  by 

fm  =  mj{\  +  ™)-  (43) 

For  example,  with  gp  —  0.5,  when  m  =  1 ,  =  0.5, 

h/hc  =  2.641;  when  m  =  9,  /m  =  0.9,  h/hc  =  7.149. 
Typically  0.5  <  m  <  2,  i.e.  0.333  </m  <  0.667,  the 
thickness  ratio  of  the  FGM  coating  is  thus 
1.882  <  h/hc  <  3.734.  It  is  also  seen  from  Fig.  15  that 
the  effect  of  different  percolation  limits  gp  on  h/fu  is 
very  small.  Consequently,  the  gp  =  0.5  curve  provides 
a  good  approximation  for  predicting  the  thickness  of 
the  FGM  coating. 

6.  CONCLUDING  REMARKS 

Functionally  graded  coatings  are  potentially  very 
attractive  to  a  number  of  applications  including  as 
thermal  barrier  coatings  and  wear  resistant  coatings. 
Compared  with  pure  ceramic  coatings,  FGM 
coatings  can  have  relatively  high  hardness  and 
oxidation  resistance  at  the  surface,  but  with  much 
lower  thermal  residual  stresses.  The  material  grada¬ 
tion  in  such  coatings  provides  a  new  dimension  in 
optimizing  the  coating  performance.  However,  due  to 
material  gradation,  the  FGM  coating  is  inherently 
heterogeneous,  and  the  microstructure-performance 
relationship  is  quite  complex  since  more  material 
parameters  come  into  play.  It  is  therefore  critical  to 
develop  mechanics  models  that  can  guide  both  the 
microstructural  design  of  a  functionally  graded 
material  and  the  design  of  components  made  of 
FGM. 

In  order  to  realize  the  potential  of  FGM  coatings, 
this  paper  focuses  on  delamination  cracking  and 
steady-state  heat  transfer  in  the  FGM  coating/metal 
substrate  system,  aiming  to  gain  a  better  understand¬ 
ing  of  the  microstructure-performance  relationship. 
Both  the  coating  and  the  substrate  are  taken  as 
linear-elastic,  with  the  coating  gradation  character¬ 
ized  by  position-dependent  Young’s  modulus  and 
thermal  expansion  coefficient.  An  edge  delamination 
crack  in  the  coating  is  analyzed  first  using  a  simple 
beam  theory,  with  the  crack  plane  parallel  to,  but  not 
necessarily  on,  the  coating/substrate  interface.  For¬ 
mulae  for  the  energy  release  rate  and  stress  intensity 
factors  of  the  crack  are  obtained  for  a  FGM  coating 
with  arbitrary  gradation;  their  dependence  on  the 
thickness  ratio  hjH,  crack  location  hi/h  and 
gradation  exponent  n  is  illustrated  for  FGM  coating 
with  Young’s  modulus  and  thermal  expansion 


coefficient  obeying  a  power-law  type  position 
dependence. 

Buckle-driven  delamination  in  a  FGM  coating  on 
an  infinitely  deep  metal  substrate  is  analyzed  based 
on  a  plate  theory;  the  critical  temperature  drop  at  the 
onset  of  buckling  and  the  energy  release  rate  of  the 
delamination  crack  are  determined  in  terms  of  the 
position-dependent  thermomechanical  properties  of 
the  coating  and  the  location  of  the  crack.  It  is  found 
that  the  use  of  FGM  coating  can  significantly  reduce 
the  likelihood  of  local  buckling  and  the  fracture 
driving  force,  since  the  thermal  residual  compressive 
stresses  in  a  FGM  coating  can  be  much  lower. 
However,  by  performing  a  steady-state  heat  transfer 
analysis,  it  is  found  that  when  the  FGM  coating  is 
used  as  a  thermal  barrier  coating,  the  thickness  of  the 
coating  needs  to  be  increased  in  order  to  provide  the 
same  degree  of  protection  to  the  metal  substrate. 

The  delamination  cracking  analysis  conducted  in 
this  study  is  based  on  a  number  of  assumptions.  For 
example,  the  crack  plane  is  taken  to  be  parallel  to  the 
interface.  In  reality,  however,  the  crack  may  be  tilted. 
Further,  plasticity  of  metal  in  the  FGM  coating  may 
play  a  very  important  role  in  delamination  cracking; 
such  an  effect  is  neglected  in  the  present  study  but  will 
be  quantified  in  future  studies.  Although  under 
far-field  uniform  tensile  straining  the  crack  driving 
force  is  found  to  be  larger  when  the  crack  is  at  the 
interface,  the  true  location  of  the  crack  is  likely  to  be 
away  from  the  interface  since  the  toughness  of  the 
FGM  coating  is  also  higher  near  the  interface. 
Modeling  the  position-dependent  fracture  toughness 
of  a  FGM  coating  is  thus  crucial. 
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CRACK  BRIDGING  IN  FUNCTIONALLY  GRADED 
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Abstract — A  crack  bridging  analysis  is  carried  out  to  predict  crack  propagation  in  coatings  made 
of  functionally  graded  materials  (FGM).  The  FGM  coating  is  taken  to  be  a  ceramic/metal  composite 
with  its  gradation  characterized  by  the  local  volume  fractions  of  metal  and  ceramic  phases.  Fracture 
in  the  FGM  coating  is  resisted  by  the  plastic  deformation  of  metal  ligaments  in  the  crack  wake  that 
bridge  the  crack ;  the  crack  bridging,  however,  is  not  uniform.  A  position-dependent  crack  bridging 
model  is  developed  taking  into  account  the  coating  gradation  and  metal  plasticity.  The  model  is 
subsequently  used  in  a  finite  element  analysis  to  predict  the  reduced  fracture  driving  force.  It  is 
found  that  crack  bridging  in  the  FGM  coating  can  significantly  reduce  the  crack  tip  stress  intensity. 
It  is  also  found  that  coating  gradation  has  a  strong  influence  on  the  fracture  driving  force  and  the 
crack  length  at  arrest.  The  present  finite  element  model  can  be  extended  readily  to  study  the  effect 
of  large-scale  plastic  deformation  on  crack  growth  in  a  FGM  coating.  ©  1997  Elsevier  Science  Ltd 


1.  INTRODUCTION 

Functionally  graded  materials  (FGMs)  have  the  potential  to  enjoy  a  wide  range  of  thermal 
and  structural  applications,  including  thermal  gradient  structures,  wear  and  corrosion 
resistant  coatings  and  metal/ceramic  joining  (Takahashi  and  Hashida,  1990;  Mortensen 
and  Suresh,  1995).  To  establish  the  fundamental  relationship  between  material  gradation 
and  thermomechanical  properties  of  FGMs,  extensive  studies  have  been  carried  out  of  the 
effective  properties  (Aboudi  et  al.  1994;  Dvorak  and  Zuiker,  1994;  Aboudi  et  al.  1995), 
deformation  and  stress  distributions  (Williamson  et  al.,  1993;  Drake  et  al.,  1993;  Gian- 
nakopoulos  et  al,  1995;  Finot  and  Suresh,  1996)  and  fracture  in  functionally  graded 
materials  (Erdogan,  1985;  Delale  and  Erdogan,  1983,  1988;  Erdogan  and  Wu,  1993; 
Chen  and  Erdogan,  1996).  In  particular,  the  fracture  mechanics  analysis  of  Erdogan  and 
associates  has  been  extended  by  Noda  and  Jin  to  include  thermal  load  (Jin  and  Noda, 
1993;  Noda  and  Jin,  1993;  Jin  and  Noda,  1994),  and  by  Bao  and  Wang  (1995)  and 
Bao  and  Cai  (1996)  to  include  the  effects  of  different  material  gradation.  Experimental 
observations  of  cracking  in  FGMs  have  been  made  by  several  investigators,  including 
surface  crack  initiation  in  thermal  barrier  coatings  (Kokini  and  Takeuchi,  1994)  and 
multiple  cracking  in  a  NiAl-Al203  FGM  layer  under  bending  (Lannutti,  1994). 

Ceramic,  metal  and  polymer  coatings  have  been  widely  used  in  various  industries, 
from  aerospace  to  automobile  to  microelectronics,  to  enhance  the  performance.  In  particu¬ 
lar,  ceramic  coatings  have  been  developed  for  use  in  thermal  gradient  structures,  automotive 
engines  and  cutting  and  grinding  tools  to  protect  the  surfaces  from  melting,  wear,  corrosion 
and  oxidation.  However,  due  to  thermal  expansion  mismatch  between  the  coating  and 
substrate,  multiple  cracking  and  coating  spalling  can  occur  upon  experiencing  thermal 
cycling  and  mechanical  loading,  causing  premature  failure  of  the  component.  The  use  of 
functionally  graded  coatings  has  the  potential  to  simultaneously  reduce  thermal  expansion 
mismatch,  increase  interface  bonding  strength,  and  enhance  coating  toughness.  To  realize 
this  potential,  a  systematic  micromechanics  study  was  carried  out  by  Bao  and  Wang 
(1995)  of  multiple  cracking  in  functionally  graded  ceramic/metal  coatings.  Systematic  finite 
element  calculations  were  made  for  the  energy  release  rate  of  the  cracks  in  the  coating  as 

f  Correspondence  address:  Department  of  Mechanical  Engineering,  The  Johns  Hopkins  University,  125 
Latrobe  Hall  3400  N.  Charles  St.,  Baltimore  MD  21218-2686,  U.S.A. 
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Fig.  1.  A  schematic  of  bridged  multiple  cracks  in  a  functionally  graded  coating/metal-substrate 

system. 


determined  by  the  coating  gradation,  crack  length,  and  the  crack  density ;  both  mechanical 
and  thermal  loads  were  considered.  However,  in  their  study,  the  influence  of  crack  bridging 
and  metal  plasticity  in  the  functionally  graded  coatings  was  neglected. 

Crack  bridging  is  a  very  common  phenomenon  in  the  fracture  of  composites  (Bao  and 
Suo,  1992).  For  example,  in  a  unidirectional  ceramic  fiber  reinforced  metal  matrix 
composite,  fibers  can  bridge  the  matrix  crack,  reducing  the  crack-tip  stress  intensity,  and 
increasing  significantly  the  fatigue  life  of  the  component  (Walls  et  al.,  1993;  Bao  and 
McMeeking,  1994,  1995).  When  a  ceramic/metal  FGM  coating  is  of  the  inclusion/matrix 
type,  the  metal  particles  in  the  ceramic  matrix  can  bridge  the  crack,  thus  enhance  the 
fracture  resistance  of  the  coating.  The  fracture  behavior  of  ductile-particle  reinforced  brittle 
matrix  composites  has  been  analyzed  (e.g.,  Bao  and  Hui,  1990;  Bao  and  Zok,  1993)  based 
on  uniform  crack  bridging  models.  However,  crack  bridging  in  a  functionally  graded 
coating  is  non-uniform,  since  the  local  volume  fraction  of  the  metal  phase  changes  with 
position  in  the  coating. 

In  this  study,  crack  bridging  in  the  FGM  coating  is  analyzed  using  a  position-dependent 
crack  bridging  law.  Considered  is  multiple  cracking  in  a  ceramic/metal  FGM  coating 
perfectly  bonded  to  a  homogeneous  metal  substrate,  as  depicted  in  Fig.  1.  The  coating  is 
taken  to  be  a  composite  with  the  local  volume  fraction  of  metal  varying  through  the  coating 
thickness  according  to  a  power-law  type  relation.  Both  the  FGM  coating  and  the  metal 
substrate  are  taken  as  linear  elastic ;  plastic  deformation  of  the  metal  phase  in  the  coating 
is  included  in  the  model  through  the  crack  bridging  analysis.  Systematic  finite  element 
predictions  are  made  for  the  fracture  driving  force  of  bridged  multiple  cracks  in  a  FGM 
coating  as  determined  by  the  crack  bridging  characteristics,  coating  gradation,  the  crack 
length  and  spacing,  and  the  applied  load.  Assuming  the  intrinsic  toughness  of  the  FGM 
coating  is  that  of  the  ceramic  phase,  the  onset  of  unstable  crack  growth  and  the  crack  arrest 
are  also  analyzed. 


2.  COATING  GRADATION  AND  THE  NON-UNIFORM  CRACK  BRIDGING 

Based  on  the  performance  requirements  for  thermal  barrier  coatings  and  coatings  in 
wear-related  applications,  in  this  study,  only  coatings  that  are  ceramic-rich  near  the  surface 
and  metal-rich  near  the  coating/substrate  interface  are  considered.  For  convenience,  the  y- 
axis  is  set  along  the  coating  thickness  direction,  the  x-axis  lies  within  the  coating/substrate 
interface,  as  illustrated  in  Fig.  1.  At  any  position  y  in  the  ceramic/metal  coating,  the  local 
volume  fraction  of  ceramic  is  assumed  to  obey  a  pure  power-law  relation  g(y) 

g(y)  =ly/hy  (!) 

where  the  exponent  n  is  a  material  parameter.  The  gradation  given  in  (1)  implies  that  the 
coating  always  has  100%  ceramic  at  the  coating  surface  (i.e.,  g{h)  =  1)  and  100%  metal  at 
the  interface  (i.e.,  g(Q)  =  0).  The  local  volume  fraction  of  metal /(y)  is  given  by 
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f(y)  =  )-g(y)  =  i-iy/hT 

(2a) 

and  the  total  volume  fraction  of  metal  fm  in  the  coating  is  related  to  the  exponent  n  in  (1) 
by 


fm  =  «/(*+!)• 


(2b) 


Due  to  the  heterogeneous  nature  of  a  functionally  graded  material,  it  is  quite  difficult 
to  obtain  exact  solutions  of  the  effective  elastic  moduli  of  a  FGM  coating  in  terms  of  the 
relative  volume  fractions  and  elastic  properties  of  the  metal  and  ceramic  phases.  For  a 
FGM  coating  of  the  inclusion/matrix  type,  the  microgeometry  of  each  phase  is  usually 
irregular,  the  spatial  distribution  of  the  phases  is  not  uniform  along  the  ^-direction,  and 
the  homogenization  can  only  be  made  within  a  plane  parallel  to  the  x-axis.  Consequently, 
the  existing  micromechanical  cell  models  are  no  longer  valid  since  no  unit  cell  containing  a 
single  inclusion  can  be  used  to  represent  the  whole  composite  body.  Even  if  one  chooses  a 
unit  cell  containing  many  particles  aligned  along  the  coating  thickness  direction,  due  to  the 
finite  thickness  of  the  coating,  there  is  a  particle-size  dependence  of  the  unit  cell.  A  robust 
micromechanical  model  for  the  thermomechanical  properties  of  a  functionally  graded 
coating  is  yet  to  be  developed  (Markworth  et  al.,  1995). 

As  a  first-order  approximation,  in  this  study,  Young’s  modulus  E(y )  and  Poisson’s 
ratio  v(y)  of  the  FGM  coating  are  obtained  using  micromechanics  models  developed  for 
composites  with  homogeneously  distributed  spherical  reinforcements.  The  effective  bulk 
modulus  k(y)  of  the  coating  can  be  expressed  as  (Hashin,  1 962) 


Ky)  =  K  + 


_ g(y)(K-km) _ 

i  +  [i  -9(y)Wc-km)Kkm+lnm)] 


(3) 


where  km  and  kc  are  bulk  moduli  for  metal  and  ceramic,  respectively,  and  /r,„  is  the  shear 
modulus  for  the  metal  phase.  The  effective  shear  modulus  Ky)  of  the  coating  can  be  solved 
from  (Christensen  and  Lo,  1979) 


Ky) 


+b 


Ky)' 


+  c  =  o 


(4) 


where  A,  B  and  C  are  functions  of  g(y),  gc/gm,  vc  and  v„,  (Poisson’s  ratios  for  ceramic  and 
metal,  respectively)  given  in  Christensen  and  Lo  (1979),  and  gc  is  the  shear  modulus  for  the 
ceramic  phase.  The  effective  Young’s  modulus  E(y )  and  Poisson’s  ratio  v(y)  are  given  in 
terms  of  k(y)  and  Ky)  by 


E{y)  = 


9  k(y)Ky) 

3  k(y)+Kyr 


vOO  = 


E(y) 
2  Ky) 


-l. 


For  a  FGM  coating  that  is  ceramic-rich  near  the  surface,  cracks  perpendicular  to  the 
interface  can  initiate  at  the  coating  surface  and  propagate  into  the  coating  under  applied 
thermal  or  mechanical  tensile  loads  (Lannutti,  1994).  These  cracks  may  be  bridged  by  the 
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Fig.  2.  A  schematic  of  position-dependent  crack  bridging  in  functionally  graded  materials. 


metal  particle  in  the  FGM  coating,  as  shown  schematically  in  Fig.  2.  As  a  consequence  of 
metal  volume  fraction  variation  in  the  FGM  coating,  the  bridging  traction  on  the  crack 
surfaces  provided  by  metal  particles  is  position-dependent  (Fig.  2).  At  a  fixed  location,  the 
bridging-stress  versus  crack-opening  curve  for  a  single  ductile  particle  bridging  the  matrix 
crack  in  its  wake  is  typically  as  the  solid  line  shown  in  Fig.  3a  (Ashby  et  al.,  1989).  To 
simplify  the  analysis  and  at  the  same  time  capture  the  essence  of  the  bridging  mechanism, 
in  this  study,  an  idealized  rectilinear  bridging  law  (Bao  and  Suo,  1992)  shown  as  the  dashed 
line  in  Fig.  3a  will  be  used.  Specifically,  the  rectilinear  bridging  law  o(5)  is  characterized  by 
two  parameters  cr0  and  <50 


a  =  <j0,  0<  <5  <c>o;  (6a) 

a  =  0,  a  >  dQ.  (6b) 

In  eqn  (6),  S0  can  be  chosen  as  the  limiting  separation  of  the  crack  surfaces  beyond  which 
the  metal  particle  will  rupture,  and  the  bridging  strength  a0  can  be  given  from  the  following 
energy-balance  equation 


cr0^o  — 


(6c) 


where  the  function  <r(<5)  is  the  measured  crack  bridging  curve  shown  as  the  solid  line  in  Fig. 
3a.  In  general,  a0  depends  on  the  plastic  flow  behavior  and  shape  of  the  metal  particles, 
and  the  particle/matrix  interface  strength  (Mataga,  1989). 

The  crack  bridging  law  given  in  (6)  is  for  a  single  metal  particle  that  bridges  the  crack. 
The  crack  bridging  law  shown  in  Fig.  3b  for  the  entire  bridging  zone  is  obtained  by  smearing 
out  the  stresses  in  the  bridging  particles  according  to  the  local  volume  fraction  of  metal  in 
the  FGM  coating 


Fig.  3.  The  crack 


(a)  (b) 

bridging  law  (a)  for  a  single  metal  particle  in  the  bridging  zone ;  (b)  for  the  FGM 
coating  where /is  the  local  volume  fraction  of  the  metal  phase. 
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a  =  ab,  0  <  §  <  50 ;  (7a) 

<r  =  0,  S  >  S0  (7b) 

where  <rA(y)  =  f(y)a0  is  the  position-dependent  bridging  strength 

0*0)  =  [1  -g(y)]<To  =  [1  -(y/h)”]a0.  (7c) 


In  deriving  eqn  (7),  it  is  assumed  that  the  local  area  fraction  of  the  bridged  particles  can  be 
taken  as  the  local  volume  fraction  of  metal  phase.  While  not  true  for  all  cases,  this 
assumption  is  valid  for  a  large  class  of  composites. 


3.  MULTIPLE  CRACKING  UNDER  UNIFORM  STRAINING 

Analyzed  in  this  section  is  multiple  cracking  in  the  FGM  coating/metal  substrate 
system  owing  to  a  remotely  applied  uniform  strain  e0  in  the  x-direction,  as  illustrated  in 
Fig.  1.  For  simplicity,  the  cracks  are  taken  to  be  mode  I,  plane  strain,  parallel  and  equally 
spaced  with  spacing  2 L,  and  have  the  same  length.  For  an  uncracked  coating/substrate 
system  subjected  to  e0,  the  stress  along  the  x-direction  in  the  coating  is  given  by 


<A(k)  = 


£o-g(j) 

i-v2(y) 


(B) 


and  that  in  the  substrate  is  uniform 


= 


E0^m 

T-72 


(9) 


The  average  stress  in  the  x-direction  corresponding  to  the  applied  strain  s0  is  obtained 
readily  as 

a~  h  +  H 


HE,,, 
1  -v2 


■  + 


o  1  —  v2  (  y)  J 


(10) 


where  H  is  the  thickness  of  the  substrate.  In  what  follows,  this  average  stress  a,  often 
referred  to  as  the  “applied  stress”,  is  used  to  represent  the  effect  of  the  applied  strain  e0. 
The  energy  release  rate  ^  at  each  tip  of  the  unbridged  multiple  cracks  can  be  normalized 
to  give 


9Em 

a2h 


=  * 


a  L  50  <7 o  <7  \ 

h’h’~h’Tm’Tm’n) 


(11) 


where  i/<  is  a  non-dimensional  function  and  the  dependence  of  if/  on  EJEm,  v,  and  vm  is  left 
implicit.  Systematic  finite  element  predictions  are  made  for  the  energy  release  rate  -f/  at  the 
tip  of  the  bridged  crack  varying  a/h,  L/h,  a0/Em,  a/Em  and  n.  A  detailed  discussion  of  these 
predictions  is  given  in  Section  3.3. 

3.1.  The  finite  element  model 

Finite  element  calculations  were  carried  out  using  the  commercial  code  ABAQUS  to 
obtain  the  energy  release  rate  ^  for  bridged  multiple  cracks  in  a  FGM  coating.  The 
coating/substrate  system  depicted  in  Fig.  1  is  subjected  to  a  remotely  applied  uniform  strain 
e0  in  the  x-direction.  Owing  to  translational  symmetry  of  the  parallel  cracks,  only  a  unit 
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cell  with  width  L  containing  a  single  crack  as  shown  in  Fig.  4  needs  to  be  considered  in  the 
finite  element  calculation.  The  thickness  of  the  substrate  H  is  taken  to  be  300 h  where  h  is 
the  coating  thickness.  The  FGM  coating  is  divided  into  100  layers  with  one  element  per 
layer  in  the  thickness  direction,  each  layer  is  homogeneous,  isotropic  with  elastic  moduli 
E(yi)  and  v (>-,■)  where  y,  is  the  position  of  the  middle  plane  of  the  layer  measured  from  the 
coating/substrate  interface  (y  =  0).  Clearly,  in  this  finite  element  model  the  gradation  is 
not  continuous,  nor  is  it  in  the  real  material.  This  model  mimics  the  microstructure  of  a 
fine  multilayered  FGM  coating,  or  an  inclusion/matrix  type  FGM  coating  with  gradation 

g(y)- 

As  can  be  seen  from  Fig.  4,  the  mode  I,  plane  strain  crack  lies  on  x  =  0,  h-a  y  ^  h. 
For  —  H  <  y  <  h  —  a,  the  displacement  in  the  x-direction  ux  =  0.  The  unit  cell  is  constrained 
in  such  a  way  that  the  edge  x  =  L  has  a  uniform  displacement  in  the  x-di  recti  on  ux  =  e0L, 
where  s0  is  the  applied  strain.  The  energy  release  rate  is  calculated  with  the  crack  tip  located 
in  the  middle  of  a  layer,  as  depicted  in  Fig.  4.  Complications  due  to  the  singular  behavior 
at  the  interfaces  of  two  adjacent  layers  can  thus  be  avoided.  The  position-dependent  crack 
bridging  is  simulated  using  discrete  nonlinear  springs  attached  to  the  nodes  on  the  crack 
surfaces  and  to  the  symmetry  plane  x  =  0.  The  behavior  of  the  nonlinear  springs  is  defined 
by  the  relationship  between  the  force  and  the  relative  displacement  in  the  spring,  i.e., 


P  = 


Abab  =  Aba0[\-(yilh)n], 

0, 


when  0  <  3  <  <50 
when  3  ^  <50 


(12) 


where  Ab  is  the  element  area  (or  length  in  the  case  of  a  plane  strain  problem)  associated 
with  the  attached  spring  and  y,  is  the  location  of  the  middle  plane  of  the  layer  containing 
the  nonlinear  spring.  In  all  the  finite  element  calculations,  fixed  values  of  Ec/Em  =  3.0, 
vc  =  0.2,  vm  =  0.3,  30/h  =  0.05  were  used.  An  accuracy  check  of  this  finite  element  crack 
bridging  model  has  been  carried  out  as  discussed  below. 

3.2.  Validation  of  the  finite  element  crack  bridging  model 

In  order  to  check  the  accuracy  of  the  finite  element  crack  bridging  model  discussed 
above,  a  special  case  is  considered  in  which  the  coating  is  taken  to  be  homogeneous  isotropic 
with  elastic  moduli  identical  to  that  of  the  substrate.  The  crack  bridging  law  is  taken  to  be 
rectilinear  as  defined  in  (7a)  and  (7b)  with  the  position-dependent  bridging  strength  <rb 
given  by 
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1 1  r 


(a) 


,aj 


l.3-0.3(b/a)' 


5/4 


(b) 


Fig.  5.  The  stress  intensity  factor  for  an  edge  crack  in  a  semi-infinite  homogeneous  isotropic  body 
under  (a)  remotely  applied  uniform  stress  and  (b)  a  pair  of  point  loads  on  the  crack  surfaces. 


Ob  =  (t/h)<J0 


(13) 


where  t  is  the  position  in  the  coating  measured  from  the  coating  surface,  and  h  is  the  coating 
thickness.  Evidently,  eqn  (13)  is  a  special  case  of  (7c)  with  t  =  h—y  and  n  =  1. 

For  a  single  unbridged  edge  crack  with  length  a  in  a  semi-infinite,  homogeneous, 
isotropic  body  under  remotely  applied  tensile  stress  a  as  depicted  in  Fig.  5a,  the  stress 
intensity  factor  Ka  is  given  by  (Tada  et  al.,  1985) 

Ka  =  \M\bayfmi.  (14) 


On  the  other  hand,  the  stress  intensity  factor  of  the  crack  under  a  pair  of  point  loads  P 
applied  on  the  crack  surfaces  a  distance  b  away  from  the  coating  surface  (Fig.  5b)  is  given 
by  (Tada  et  al.,  1985) 


where 

/A  _  1.3  — 0.3(fr/fl)5/4 

W  J\-(b/a)2 


(15a) 


(15b) 


The  crack-tip  stress  intensity  Ktip  for  a  fully  bridged  crack  is  thus 


Substituting  (13)  into  (16)  and  performing  the  integration,  we  have 


^1.1215  — 


2.15  cr0  a 

7i  ah 


(16) 


(17) 


In  a  nondimensional  form,  the  energy  release  rate  <§  at  the  crack-tip  can  be  expressed  as 
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Crack  Length,  a  /  h 

Fig.  6.  Comparison  of  the  normalized  crack-tip  energy  release  rate  as  a  function  of  the  crack  length 
obtained  using  the  closed-form  solution  (18)  and  the  finite  element  crack  bridging  model. 


VEm 

a2h 


=  n 


^1.1215  — 


2.15  cr0  aY  a 
n  ahjh 


(18) 


where  Em  =  EJ{\  —  v}n)  is  the  plane  strain  Young’s  modulus  of  the  substrate. 

Plotted  in  Fig.  6  is  the  comparison  of  the  normalized  energy  release  rate  as  a  function 
of  the  crack  length  obtained  from  the  closed-form  solution  (18)  and  from  the  finite  element 
calculations  with  a0/Em  =  0.001  under  applied  stress  a/Em  =  0.001.  In  the  finite  element 
calculations,  L/h  =  100,  Hfh  =  300  is  used.  The  limiting-separation  <50  assumed  in  the  crack 
bridging  law  is  sufficiently  large  (50/h  =  0.05)  so  that  the  crack  remains  fully  bridged  under 
the  applied  load.  As  can  be  seen  from  the  comparison  shown  in  Fig.  6,  the  finite  element 
model  is  quite  accurate  in  obtaining  the  energy  release  rate  of  a  bridged  crack. 

3.3.  Finite  element  predictions 

Plotted  in  Fig.  7  is  the  normalized  energy  release  rate  yEJo2h  as  a  function  of 
the  normalized  crack  length  afh  with  n  =  1,  o0/Em  =  0.002,  L/h  =  2  under  applied  stress 
a/Em  =  0.0004.  Corresponding  results  for  multiple  cracks  in  a  pure  ceramic  coating  and 
unbridged  cracks  in  the  FGM  coating  are  also  shown  for  comparison.  It  is  clear  that  FGM 
coatings  can  significantly  reduce  the  crack  driving  force  and  crack  bridging  can  make  the 
driving  force  reduction  even  more  pronounced.  In  particular,  for  the  bridged  cracks,  the 
energy  release  rate  versus  crack  length  curve  shows  a  peak  indicating  that  the  cracks  may 
arrest  at  a  certain  crack  length  when  the  fracture  driving  force  falls  below  the  intrinsic 
toughness  of  the  coating  and  thereafter,  an  increased  applied  load  is  required  to  propagate 
the  crack  further  (stable  crack  growth).  This  scenario  will  be  discussed  in  more  detail  in 
Section  4.  It  is  also  noteworthy  that  for  the  bridged  crack,  the  energy  release  rate  ^  at  the 
crack  tip  goes  to  zero  at  about  a  =  0.6 h.  This  is  due  to  the  fact  that  the  applied  load  is 


Crack  Length,  a/h 

Fig.  7.  Comparison  of  the  normalized  energy  release  rates  as  a  function  of  the  crack  length  for  (i) 
unbridged  cracks  in  a  pure  ceramic  coating ;  (ii)  unbridged  cracks  in  a  FGM  coating ;  and  (iii) 

bridged  cracks  in  a  FGM  coating. 


relatively  low,  and  the  crack  bridging  is  quite  strong ;  the  crack  is  always  fully  bridged  since 
the  assumed  limiting-separation  <50  is  large. 

To  uncover  the  influence  of  coating  gradation  on  the  fracture  driving  force,  systematic 
predictions  are  made  for  the  energy  release  rate  corresponding  to  different  coating 
gradation  g(y).  Shown  in  Fig.  8  are  ^EjcEh  vs  a/h  curves  for  L/h  =  2,  cr0 /Em  =  0.002, 
a/Em  =  0.0005  for  n  =  0.5,  1  and  2.  As  expected,  there  is  a  strong  dependence  of  on  the 
value  of  gradation  exponent  n :  a  large  value  of  n  leads  to  a  low  energy  release  rate  <$.  Note 
that  a  large  n  implies  a  high  metal  content  in  the  coating,  as  indicated  by  eqn  (2b). 

The  effect  of  bridging  strength  <t0  is  illustrated  in  Fig.  9  in  which  gEm/<j2h  is  plotted  as 
a  function  of  a/h  for  n  =  1,  L/h  =  2,  a/Em  =  0.0004  for  a0/Em  =  0.001,  0.002,  0.003.  The 
trend  is  clear :  a  high  bridging  strength  a0  gives  a  low  fracture  driving  force.  In  fact,  with  a 
fixed  crack  opening  6  at  the  coating  surface  the  reduction  in  crack-tip  energy  release  rate 
due  to  plastic  deformation  of  the  bridging  metal  particles  is  simply  proportional  to  a0.  Note 
also  that  the  crack  length  at  which  the  energy  release  rate  rS  approaches  zero  decreases  with 
increasing  <r0- 

Presented  in  Fig.  10  is  <$Em/o2h  as  a  function  of  the  crack  length  a/h  for  n  =  1, 
a/Em  =  0.0004,  Co /Em  =  0.002  for  L/h  =  2,  5,  10  and  50.  The  energy  release  rate  increases 
with  crack  spacing,  or  in  other  words,  the  fracture  driving  force  for  each  crack  is  reduced 
as  more  and  more  cracks  developed  (i.e.,  L/h  becomes  smaller).  It  is  interesting  to  note 
that  when  the  crack  length  is  small  (say,  a/h  <  0.2),  the  energy  release  rate  is  essentially 
independent  of  crack  spacing.  It  is  also  clear  that  when  crack  spacing  L/h  becomes  large 
(say,  L/h  >  50),  the  interaction  between  the  parallel  cracks  diminishes,  and  each  crack 
behaves  like  an  isolated  crack.  That  is  the  reason  behind  the  saturating  trend  of  the  effect 
of  L/h  on  energy  release  rate  as  shown  in  Fig.  10. 


Normalized  Energy  Release  Rate,  ?Em/CT2h  Normalized  Energy  Release  Rate,gEm/o 
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Crack  Length,  a/h 

Fig.  8.  The  normalized  energy  release  rate  as  a  function  of  the  crack  length  for  FGM  coatings  with 
n  =  0.5,  1.0  and  2.0,  Ljh  =  2,  aJEm  =  0.002  under  applied  stress  a/Em  =  0.0004. 


Crack  Length,  a/h 

Fig.  9.  The  normalized  energy  release  rate  as  a  function  of  the  crack  length  for  FGM  coatings  with 
n  =  1,  L/h  =  2,  o-0 =  0.001,  0.002  and  0.003  under  applied  stress  <r/Ew  =  0.0004. 
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Crack  Length,  a/h 

Fig.  10.  The  normalized  energy  release  rate  as  a  function  of  the  crack  length  for  a  FGM  coating 
with  n  =  1,  L/h  =  2,  5,  10  and  50,  a0IE„,  =  0.002  under  applied  stress  o/Em  =  0.0004. 


Figure  1 1  shows  the  normalized  energy  release  rate  3 EJa2h  as  a  function  of  the  crack 
length  a/h  for  n  =  1,  Llh  =  2,  a0/Em  =  0.002  for  different  applied  load  a/Em  =  0.0003, 
0.0004, 0.0005.  Evidently,  the  fracture  driving  force  3  increases  with  applied  load.  However, 
the  crack-tip  energy  release  rate  3  does  not  scale  with  tr2  as  for  an  unbridged  crack  due  to 
the  nonlinear  feature  of  the  crack  bridging  and  the  interplay  between  crack  bridging  and 
the  applied  load.  Otherwise,  all  the  curves  shown  in  Fig.  1 1  should  collapse  onto  one  curve. 

The  yEjcfh  vs  a/h  curves  shown  in  Figs  8-1 1  indicate  that,  under  the  given  applied 
far-field  uniform  strain  e0,  the  energy  release  rate  3  at  the  crack-tip  may  become  zero  at 
certain  crack  length.  This  can  be  understood  from  the  following  energy  equation 

^  =  ^app.-^  (19) 

where  ^appl  and  3b  are  the  energy  release  rates  due  to  the  applied  strain  e0  and  crack 
bridging,  respectively.  Owing  to  the  functional  gradation  of  the  coating,  for  ajh  >  0.2,  {#appl 
does  not  change  much  with  crack  length  as  can  be  seen  from  Fig.  7,  but  3b  increases  as  the 
crack  length  a  increases.  At  certain  crack  length  say,  a  —  ac,  3b  may  become  the  same  as 
^appb  which  implies  that  at  a,  the  crack-tip  energy  release  rate  3  is  zero. 

4.  STABLE  CRACK  GROWTH  IN  A  FGM  COATING 

As  mentioned  earlier,  the  fracture  driving  force  vs  crack  length  curve  for  bridged 
multiple  cracks  in  a  functionally  graded  coating  has  a  peak  indicating  that  a  crack  may 
initiate  at  certain  crack  length  and  arrest  at  another  under  a  given  applied  load.  Based  on 
the  finite  element  results  for  the  crack-tip  energy  release  rate  discussed  in  the  last  section, 
and  the  intrinsic  toughness  of  the  coating,  the  crack  length  at  the  onset  of  unstable  crack 
growth  and  the  crack  length  at  arrest  can  be  quantified. 
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Crack  Length,  a/h 

Fig.  11.  The  normalized  energy  release  rate  as  a  function  of  the  crack  length  for  a  FGM  coating 
with  n  =  1,  L/h  =  2,  <r0/Em  =  0.002  under  applied  stresses  a/E„,  =  0.0003,  0.0004  and  0.0005. 


4. 1 .  Intrinsic  toughness  of  a  FGM  coating 

As  a  ceramic/metal  composite,  the  toughness  of  a  FGM  coating  is  in  general  determined 
by  the  toughness  of  both  the  ceramic  and  the  metal  phases.  However,  compared  with  the 
metal  phase,  the  ceramic  phase  in  a  FGM  coating  experiences  larger  stresses  under  a 
uniform  applied  strain  field,  but  its  toughness  is  much  lower.  Therefore,  the  ceramic  phase 
in  the  coating  always  fractures  first.  In  this  study,  we  assume  that  in  the  ceramic-rich  region 
of  the  coating,  the  cracks  are  in  the  ceramic  matrix ;  metal  particles  contribute  to  the  coating 
toughness  only  through  crack  bridging,  as  illustrated  in  Fig.  12.  Further,  in  the  metal-rich 
region,  the  crack  length  is  so  defined  such  that  the  fractured  ceramic  particles  lying  within 
the  major  crack  plane  are  included  (Fig.  12).  Specifically,  the  crack  tip  is  assumed  to  be 
located  at  the  breaks  of  the  ceramic  particles,  and  the  metal  matrix  in  the  crack  plane 
behind  the  crack-tip  is  modeled  as  the  bridging  material.  This  approach  is  analogous  to  the 
cohesive  zone  model  used  by  Rice  (1968)  for  crack  growth  in  a  metal  under  small-scale 
yielding  conditions.  We  thus  define  the  intrinsic  toughness  of  a  FGM  coating  as  that  of  the 
ceramic  phase ;  the  contribution  of  metal  phase  to  the  coating  toughness  is  reflected  by  the 
reduction  in  the  crack-tip  energy  release  rate  f#  due  to  crack  bridging.  From  the  viewpoint 
of  energy  balance  concerning  crack  growth,  both  approaches,  i.e.,  either  counting  bridging 
as  a  toughness  enhancer  or  as  a  driving  force  reducer,  are  equivalent. 

Assume  that  the  toughness  of  the  ceramic  phase  is  Tc.  For  typical  ceramic/metal 
combinations  for  advanced  composite  materials,  the  normalized  intrinsic  toughness  of  the 
coating  YJEJi  can  range  from  1.0  x  10~7  to  1.0  x  10“5  (Ashby  and  Jones,  1980).  To  gain 
insight  into  the  stable  crack  growth  in  a  FGM  coating,  in  the  following  discussions,  the 
material  system  consisting  of  SiC/Ti  FGM  coating  on  a  Ti  substrate  will  be  considered  as 
an  example.  Since  for  silicon  carbide  (SiC),  toughness  Kc  =  3  MPa^/m,  Young’s  modulus 
Ec  =  450  GPa,  and  for  titanium  (Ti),  Em=  116  GPa  (Ashby  and  Jones,  1980),  we  have 
TJEJi  =  1.66  x  10-7  for  a  SiC/Ti  coating  with  1  mm  thickness.  For  material  systems  other 
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Ceramic-Rich  Region 


Ceramic  Metal 


Metal-Rich  Region 

Fig.  12.  A  schematic  showing  crack  bridging  in  different  regions  in  the  FGM  coating.  The  crack  tip 
is  always  in  the  ceramic  phase  and  the  intrinsic  toughness  of  the  coating  is  assumed  to  be  that  of 

the  ceramic. 


than  the  SiC/Ti  considered  here,  the  same  approach  can  be  taken  in  predicting  the  stable 
crack  growth  in  a  FGM  coating. 

4.2.  Stable  crack  growth  in  a  FGM  coating 

Unstable  crack  growth  is  a  common  feature  for  brittle  materials  such  as  ceramics  of 
which  the  toughness  is  a  constant  independent  of  the  crack  length.  However,  the  fracture 
driving  force  is  usually  a  monotonic  increasing  function  of  crack  length.  Hence,  once  the 
crack  starts  to  grow  the  whole  structure  or  component  ruptures  since  the  fracture  driving 
force  is  always  higher  than  the  toughness  after  the  crack  growth  initiates  at 

9  =  rc  (20) 

where  9  is  the  fracture  driving  force  and  Tc  is  the  material  toughness. 

In  contrast  to  unstable  crack  growth,  stable  crack  growth  implies  that  an  increased 
applied  load  is  necessary  to  propagate  the  crack  further.  For  bridged  multiple  cracks  in  a 
FGM  coating,  under  a  given  applied  load,  the  fracture  driving  force  9  increases  with 
crack  length  at  the  beginning,  reaches  the  peak  value,  then  decreases  with  crack  length. 
Consequently,  even  if  the  crack  starts  to  grow  under  a  given  applied  load,  the  crack  may 
arrest  at  certain  crack  length.  Once  the  crack  arrests,  an  increase  in  the  applied  load  is 
required  for  further  crack  growth. 

Shown  in  Fig.  13a  is  a  typical  energy  release  rate  9EJa2h  vs  crack  length  ajh  curve 
along  with  the  intrinsic  toughness  of  the  coating  YcEJa2h  represented  by  the  horizontal 
line,  which  intersects  the  driving  force  curve  at  crack  lengths  a,  and  as.  These  two  points 
have  special  meaning  and  importance :  a, ■  corresponds  to  the  critical  initial  crack  length 
beyond  which  crack  growth  starts  and  as  is  the  crack  length  at  arrest.  Under  a  given  applied 
load,  if  the  initial  crack  size  is  smaller  than  a:  the  crack  will  not  grow  because  the  fracture 
driving  force  is  lower  than  the  toughness.  On  the  other  hand,  if  the  initial  crack  size  is 
between  a,  and  as,  the  crack  will  grow,  and  eventually  arrest  at  a  =  as. 


Normalized  Energy  Release  Rate,  ^Em/o2h  Normalized  Energy  Release  Rate,  £Em/a2h 
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Crack  Length,  a/h 


Crack  Length,  a/h 

Fig.  13.  Energy  release  rate  as  a  function  of  the  crack  length  together  with  the  intrinsic  toughness 
of  the  coating,  (a)  Under  a  given  applied  load,  the  crack  starts  to  grow  if  a  >  a „  and  arrests  at 
a  =  a,,  (b)  With  different  applied  loads,  the  crack  has  different  critical  initial  crack  length  a,  and 

crack  length  at  arrest  as. 
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To  further  illustrate  the  trends  in  crack  growth,  in  Fig.  13b  curves  of  energy  release 
rate  <& EJa2h  vs  crack  length  a/h  are  plotted  together  with  the  intrinsic  toughness  of  the 
coating  rc£ja2h  represented  by  the  horizontal  lines  for  n  =  1 ,  L/h  =  2,  a0/Em  =  0.002  for 
ajEm  =  0.0004,  0.0005  and  0.0006.  Note  that  the  intrinsic  toughness  of  the  coating  is 
assumed  to  be  TJEJi  =  1 .66  x  10-7  and  different  horizontal  lines  representing  TcEjo2h 
correspond  to  different  values  of  the  applied  load  used  in  calculating  the  energy  release  rate 
<&.  Evidently,  when  the  applied  load  is  low,  say,  <y/Em  =  0.0004,  the  whole  curve  of  fracture 
driving  force  ^  falls  below  the  toughness  curve  Tc;  crack  growth  is  thus  prohibited.  Under 
a  larger  applied  load,  say,  a/Em  =  0.0005,  the  multiple  cracks  grow  when  their  length 
a  ^  0.05 h,  and  arrest  at  about  a  =  0.54/2.  At  an  even  higher  load,  a/Em  =  0.0006,  the  cracks 
start  to  grow  at  a  x  0.01/r  and  stop  at  a  =  0.7/2.  As  expected,  beyond  certain  crack  length, 
further  crack  growth  in  the  coating  requires  an  increase  in  the  applied  load. 

As  discussed  in  the  previous  section,  the  fracture  driving  force  for  bridged  multiple 
cracks  in  the  FGM  coating  depends  on  coating  gradation  as  well  as  the  crack  bridging 
characteristics.  In  general,  under  a  given  applied  load,  the  crack  length  at  arrest  as  increases 
with  decreasing  n  and  o0/Em,  as  can  be  seen  from  the  curves  in  Figs  8  and  9.  The  influence 
of  n  and  a0/Em  on  the  critical  initial  crack  length  a,  is  found  to  be  relatively  insignificant. 

5.  CONCLUDING  REMARKS 

Functionally  graded  ceramic/metal  coatings  can  have  fracture  toughness  much  higher 
than  the  pure  ceramic  coating  due  to  the  plastic  deformation  of  the  metal  phase  in  the 
coating.  Thus,  cracks  initiated  at  the  coating  surface  and  propagating  perpendicular  to  the 
coating/substrate  interface  must  overcome  greater  resistance.  Specifically,  in  the  ceramic- 
rich  region  of  the  coating,  metal  ligaments  in  the  crack  wake  bridge  the  crack ;  in  the  metal- 
rich  region,  plastic  deformation  of  the  metal  matrix  near  the  crack  tip  consumes  energy. 
To  quantify  the  effect  of  crack  bridging  and  plastic  deformation  on  crack  growth  in 
the  FGM  coating,  in  this  article,  a  finite  element  model  for  multiple  cracking  in  FGM 
coating/metal  substrate  systems  is  developed  based  on  a  position-dependent  crack  bridging 
law.  The  influence  of  material  and  mechanics  parameters  such  as  coating  gradation,  bridging 
strength,  crack  length  and  spacing,  and  applied  load  on  the  crack  tip  energy  release  rate  is 
studied  systematically.  Though  relatively  simple,  the  crack  bridging  analysis  given  in  this 
paper  can  provide  a  basis  for  further  studies  of  the  fracture  toughness  of  functionally 
graded  materials. 

The  finite  element  results  reveal  that  the  coating  gradation  exponent  n  can  have 
significant  effects  on  the  fracture  driving  force  and  the  crack  length  at  arrest :  under  fixed 
applied  load,  a  larger  value  of  n  (i.e.,  higher  metal  content  in  the  coating)  gives  a  lower 
energy  release  rate  ^  and  a  smaller  crack  length  at  arrest  as.  Further,  with  the  same  coating 
gradation,  a  high  bridging  strength  <r0  leads  to  a  large  reduction  in  the  fracture  driving  force 
and  the  crack  length  at  arrest.  However,  coating  gradation,  crack  spacing,  and  bridging 
strength  have  very  limited  effects  on  the  critical  initial  crack  length  a ,. 

In  this  study  we  assume  that  the  ductility  <50  of  the  bridging  material  (the  separation 
of  the  crack  surfaces  at  particle  rupture)  and  the  intrinsic  toughness  of  the  coating  Tc  do 
not  change  with  position  in  the  coating.  However,  in  reality,  the  crack  bridging  parameters 
ab  and  <50,  and  the  intrinsic  toughness  Tc  may  all  be  position-dependent.  Usually  S0  is 
influenced  by  the  size  and  shape  of  the  metal  particles  which  may  change  as  the  local  volume 
fraction  of  metal  in  the  coating  varies.  Further,  the  intrinsic  toughness  Fc — defined  as  the 
toughness  of  the  ceramic  phase — may  depend  on  microcracking  in  the  ceramic  phase,  and 
the  thermal  residual  stresses.  The  precise  position-dependence  of  ab  may  also  be  more 
complicated  than  that  assumed  in  this  article  ( cb(y )  =  [1  —  (y/h)"]a0),  since  <r0  may  be 
different  in  the  ceramic-rich  region  where  the  metal  particles  are  surrounded  by  the  ceramic 
matrix,  and  in  the  metal-rich  region  where  the  metal  is  the  continuous  phase.  All  these 
issues  will  be  addressed  in  the  subsequent  studies  of  the  fracture  resistance  of  functionally 
graded  materials. 

The  finite  element  crack  bridging  model  developed  in  this  study  is  based  on  several 
assumptions  which  render  the  model  to  be  applicable  only  for  crack  lengths  a  <  0.7/2.  First, 
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the  limiting  separation  of  the  crack  bridging  is  taken  to  be  S0  =  0.05 h  which  is  quite  large. 
Consequently,  in  all  the  finite  element  calculations,  the  crack  remains  fully  bridged.  This 
may  have  overestimated  the  effect  of  crack  bridging  when  the  crack  length  becomes  large. 
Further,  the  plastic  deformation  of  metal  in  the  metal-rich  region  of  the  FGM  coating  is 
taken  into  account  through  the  position-dependent  crack  bridging  law  by  assuming  that 
the  crack  tip  is  located  at  the  fractured  ceramic  particles  in  the  crack  plane.  A  consequence 
of  this  approach  is  that  the  model  predictions  are  less  accurate  for  crack  lengths  much 
larger  than  0.6/z.  This  model,  similar  to  the  cohesive-zone  model  used  by  Rice  (1968),  is 
valid  only  if  the  plastic  yielding  is  of  small  scale.  For  large-scale  yielding  cases,  the  present 
finite  element  model  can  be  extended  readily  to  fully  reflect  the  effect  of  plastic  deformation 
on  crack  growth  in  a  FGM  coating  by  taking  the  metal-rich  region  of  the  FGM  coating 
and  the  metal  substrate  to  be  elastoplastic  materials.  Such  an  extension  is  left  for  future 
studies. 
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